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ABSmCT 


As part of its policy of supporting research and development programs 
which reside on the frontier of power technology, the National Aeronautics 
and Space Administration has sponsored work in gaseous fueled reactors and 
plasma research. The original thrust of the NASA sponsored research, aimed 
toward development of a space propulsion engine, led to two gas-core reactor 
concepts — the light bulb and the coaxial flow nuclear reactor concepts. 
Although budgetary and policy factors terminated the development of nuclear 
powered propulsion engines, the concept of a TJFg fueled gas core reactor was 
shown to be very attractive for several other applications. Studies have 
shown that UFg fueled reactors can be quite versatile with respect to power, 
pressure, operating temperature, and the modes of power extraction. Pos- 
sible cycles Include Brayton cycles, Eankine cycles, MHD generators, and 
thermionic diodes. Recent results of research on the pumping of lasers 
by fission fragment interactions with a laser gas mixture indicate the pos- 
sibility of power extraction in the fom of coherent light. Another poten- 
tial application of the gas core reactor is its use for nuclear waste dis- 
posal by nuclear transmutation. 

Circulating-fuel breeder reactors with continuous fission-product 
removal and simplified on-site reprocessing have a potential for generating 
cheaper power and greatly reducing the safety and safeguards problems asso- 
ciated with fuel transportation. In addition, uranium hexafluoride gas 
fueled breeder reactors may "bum up” most, or possibly all, of their long- 
lived radioactive waste products. Because of the highly favorable economics, 
greatly simplified on-site fuel processing, enhanced safety (no fuel melt- 
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down possible), and considerably reduced safeguards and radwaste problems, 
the uranium-hexafluoride breeder reactor may well be the solution to the 
problems facing the nuclear industry. 

233 

The reactor concept analyzed was a TH?g core surrounded by a molten 
7 

salt (Li F, beF 2 , ThFg) blanket. Nuclear survey calculations were carried 
out for both spherical and cylindrical geometries. The nuclear analysis 
techniques were assessed by comparison with a previous check-point calcula- 
tion. The core diameter was approximately one meter. The blanket thick- 
ness ranged from 60 to 130 cm. A breeding ratio of approximately l.ld and 
a critical mass about 379 kg were found. Further neutronic calculations 
were made to assess the effect on the critical mass, breeding ratio, and 
spectrum of substituting moderator. Be or G, for the molten salt in the 
blanket . 

The major problem of the materials selection was to find materials 
compatible with the OFg in the core region and the molten salt in the blanket 
region. The materials selected for use in the reactor system were Monel 
and modified Hastelloy-N. Modified Hastelloy-N is an alloy developed spe- 
cifically for use with molten fluoride systems. Its major constituent is 
nickel. The ductility of modified Hastelloy-N enables its use together as 
a liner between the core and blanket and also the containment structure for 
the blanket. Another problem is the selection of turbine blade materials. 

The high temperature of the UFg gas poses severe problems as far as corro- 
sion and strength are concerned, 

/ 

Thermodynamic cycle calculations were performed for the Rankine cycle. 
Calculations were also done for a cycle with no reheater to determine how 
the cycle efficiency is affected by the reheater. Because material limita- 
tions of the duct walls limited UFg temperature to less than 1660°R C921,89”IC), 
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the average temperature at the outlet of the reactor was chosen to he 1560°R 
(866,33°K). The reactor outlet pressure was set at 1450 psia (99.97 bars) 
which is approximately the pressure required from the core physics calcula- 
tions to obtain a critical reactor, 

The minimum temperature difference in the regenerative heat exchanger, 
ATg_y, was specified as 50 °R (28°K). The boiler feed pump and turbine effi- 
ciencies were taken to be 0.88 and the condenser pressure was maintained at 
21.76 psia (1.5 bars)^. The pressure loss across the reactor was assumed to 
be 14.89 psl (1,03 bars) and the pressure loss across each side of the re- 
generator and reheater was taken to be 7.445 psi (.51 bar). In these calcu- 
lations, the pressure loss through the pipes was neglected because the plant 
design was not complete and therefore the length of the pipes between com- 
ponents was unknown. 

Initially, cycle parameters were calculated without reheating. This 
gave a plant efficiency of 39.03 percent. Then reheat was added between 
stages of the turbine using the heat generated in the blanket. It was ap- 
parent that an optimum pressure existed for removing the UFg from the high 
pressure turbine and, consequently, p2 was varied. Overall plant efficien- 
cies as a function of the high pressure turbine outlet pressure were cal- 
culated, The maximum efficiency was 41.44% at a pressure of 435,12 psia 
(30 bars). The efficiency begins to increase for pressures greater than 
720 psia (50 bars). However, higher pressures cannot be used because calcu- 
lations indicate that for these higher pressures the reheat temperatures ex- 
ceed 1560°R (866®K) which was not allowed because of material limitations. 

A conceptual mechanical design is presented along with a system layout 
for a 1000 MW stationary power plant. Further work is required before a 
final design can be completed. 
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In summary, the advantages of the GCBR include : 

(1) high efficiency, 

(2) simplified on-line reprocessing, 

(3) inherent safety considerations, 

(4) high breeding ratio, 

(5) possibility of burning all or most of the long-lived nuclear waste 
actinides, and 

(6) possibility of extrapolating the technology to higher temperatures 
and MHD direct conversion. 


I. INTRODUCTION 


As part of its policy of supporting research and development programs 
which reside on the frontier of power technology, the National Aeronautics 
and Space Administration has sponsored work in gaseous fueled reactors and 
plasma research. The original thrust of the NASA sponsored research, aimed 
toward development of a space propulsion engine, led to two gas-core reactor 
concepts - the light bulb and the coaxial flow nuclear reactor concepts.^ 
Although budgetary and policy factors terminated the development of 
nuclear powered propulsion engines, the concept of a UFg fueled gas core 
reactor was shown to be very attractive for several other applications . 

NASA has continued supporting an ongoing fissioning plasma research 
program consisting of cavity reactor criticality tests, fluid mechanics 
tests, investigations of uranium optical emission spectra, radiant heat trans- 
fer power plant studies, and related theoretical work. These studies 

have shown that UFg fueled reactors can be quite versatile with respect to 
power, pressure, operating temperature, and the modes of power extraction. 

Possible power conversion systems include Bray ton cycles, Ranklne cycles, 

12 13 14 

MHD generators, and thermionic diodes. ’ * Recent results of research 

on the pumping of lasers by fission fragment interactions with a laser gas 
mixture indicate the possibility of the power extraction in the form of co- 
herent light. Another potential application of the gas core reactor is 
its use for nuclear waste disposal by nuclear transmutation (Gas Core Actinide 
Transmutation Reactor, GCATR) . ^^»12,15,17 
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One of the major advantages of UFg reactors for power generation is the 
simplified fuel reprocessing scheme which the gaseous fuel makes possible. 

Part of the UFg from the reactor can be passed through a fractional distilla- 
tion process to remove fission products with higher or lower boiling points 
than the pure UFg. Also, because of the large mass difference, fission pro- 
ducts can be separated from the UFg in a centrifuge. The long-lived radio- 
active fission products are appropriately returned to the reactor where they 
undergo intense neutron irradiation in the reactor core, or in a region in 
the blanket near the core interface. In the latter case the irradiated fis- 
sion product mixture gradually diffuses into the core to mix with the gaseous 
fuel. During irradiation, a large fraction of the long-lived fission products 
are transmuted into short-lived radioactive isotopes. In this manner, the 
long-lived radioactive Isotopes are recycled until they are burned up. Thus, 
the UF0 reactor produces considerably less radioactive waste than other fis- 
sion reactors for disposal, and possibly none of the long-lived isotopes. 

The UFg breeder reactor is inherently safe because the conventional loss- 
of-coolant accident cannot occur, radioactive fission products are continuously 
removed from the fuel stream resulting in steady-state operation with a core 
that contains only a minimum steady-state concentration of fission products, 
and the temperature coefficient or reactivity is negative which prevents acciden- 
tal power excursions. Another important safety aspect of the UFg breeder is the 
use of U-233 Instead of Pu-239 as the fissile isotope; U-233 is considerably 
less toxic. Work at Georgia Tech has demonstrated the potential of UFg fueled 
reactors for breeding and investigated several energy conversion schemes The 
nuclear analysis method used multigroup diffusion theory. A benchmark experi- 
ment carried out at Idaho Falls was analyzed to show that the nuclear analytical 


technique was sufficiently accurate. The UF^ fueled breeder reactor was con- 
sidered spherical with U-233-Fg as fuel and beryllium and thorium in the 
blanket. A breeding ratio of 1.25-1.26 was calculated for core diameters 
varying from 1 to 5m. Fuel doubling times were a few years or less. Conven- 
tional Rankine cycle and Brayton cycle efficiencies ranged from 35% to 50%, 

A further analysis of the gas core breeder reactor has been sponsored 
at Georgia Tech under NASA Research Grant 1168 (April 6, 1975 to February 6, 
1976), This final report describes the results to date. 
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XX. BACKGROUND 


N76-19554 


Critical Experiments 

The first uranium hexafluoride gas fueled reactor built and operated 
in the Soviet Union in the 1950 's,^ is shown in Fig. XX-1, The uranitim hexa- 
fluoride gas flowed through a 148 channel beryllium moderator which was 
surrounded by a 50 cm thick graphite reflector. This reactor was critical 
with 3.34 Kg (7.36 lb) of 90% enriched uranium hexafluoride gas and oper- 
ated at a power level of 1.5 Kwt. Two boron carbide control rods were in- 
serted into the core to control the reactivity. Table IX-1 described the 
characteristics of the reactor. 

Table IX-1 Characteristics of the Soviet DFg Gas Fueled 
Reactor Operated in the 1950's^ 


Core diameter, m (ft) 

1.16 

(3.81) 

Core height, m (ft) 

1.08 

(3.51) 

Core gas volume, 1 (ft^) 

213 

(7.52) 

Weight of Be in core. Kg (lb) 

1370 

(3020) 

Weight of A1 in core. Kg (lb) 

62.8 

(138) 

Graphite reflector thickness, cm (in) 

50 

(20) 

Pitch of square lattice, cm (in) 

8 

(2.15) 

Number of channels 

148 


Critical mass of UFg, Kg (lb) 
Critical mass of U-235, Kg (lb) 

3.34 

(7.36) 

2.024 

(4.46) 


This reactor demonstrated several positive aspects of uranium hexaflu- 
oride fueled reactors. Since refueling is done continuously » the reactor 
is always operating with a "clean core," and there is no need for fuel 
fabrication. The absence of poisons in the core significantly improves the 
neutron economy and breeding ratio, and reduces the fuel doubling time. 
Besides, the reactor shows large negative temperature reactivity due to the 
large expansion coefficient of gaseous UFg in the core. 
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Figure II-l UF^ Fueled Reactor Operated in the U.S.S.R. in the 1950 's 
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In recent years, interest In uranium hexafluoride reactors has turned 
toward the cavity reactor concept in which the uranium hexafluoride gas is 
contained within a large cavity surrounded fay a neutron moderator. The 
reason for pursuing the cavity reactor approach is the higher temperatures 
which may fae possible with such reactor system. However, the earlier ap- 
proach of using channels in a solid moderator such as graphite should not 
be ruled out, since in the near future the UFg reactor with somewhat lower 
temperature may be possible. 

On May 17, 1967, a cylindrical uranium hexafluoride gaseous core cavity 

2 

reactor was made critical at Idaho Falls, Idaho. The gaseous uranium 
hexafluoride was contained inside a cylindrical aluminum tank with a diam- 
eter of 1.22 meters (4 feet) and a length of 1.09 meters (3.58 feet) as 
illustrated in Fig. II-2, The surrounding moderator-reflector tank was a 
cylinder with an Inside diameter of 1.83 meters (6 feet) and a length of 
1.22 meters (4 feet). The outside diameter was 3.66 meters (12 feet) and 
its length was 3,05 meters (10 feet). The annular reflector region was 
filled with heavy T'rater to a nominal thickness of 0.89 meter (3 feet) and 
contained a 10.2 centimeter (4 inch) thick beryllium slab placed 6.7 cm 
(2,64 inches) from the cavity wall. The UFg fuel \ms 93% enriched in the 
U-235 isotope, and the experiment had a critical mass of 21.3 Kg (47 lb) of 
U-235 . 

This UFg cavity reactor critical facility used a split table arrangement 

with one of the tanks consisting of an end reflector Tdiich could be pulled 

* 

back for easy access to the cavity region. The cavity was contained \ri.thin 
the larger fixed tank and was 1,22 meters (4 feet) long and 1,83 meters (6 
feet) in diameter. The surrounding tanks contained heavy water which served 
as the reflector-moderator. The reflector region was nominally 88.9 cm (35 in) 
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thick. The reactor was controlled with stainless steel clad boron carbide 
rods which were insetted into the end of tdie fixed tank. 



Tig. II-2 First U.S. UFg Reactor Critical Experiment^ 


During the latter half of 1969 a spherical cavity reactor experiment 
was constructed and operated with enriched uranium hexafluoride gas at the 
National Reactor Test Site in Idaho. The uranium hexafluoride gas was 
contained within a tank of aluminum surrounded by a heavy water moderator 
blanket. The tank was spherical, consisting of two spun hemispheres welded 
together, with dimensions of 1.27 meters (4.18 feet) outside diameter and 
0.635 centimeter (0.25 inch) wall thickness. The 90.72 Kg (200 lb) tank 
was designed for a maximum internal pressure of 4.46 bars (64,7 psi) and a 
maximum temperature of 400“K (260‘’F). In order to keep the uranium hexa- 
fluoride gas from condensing, heated air flowed around the core between 
the outside of the spherical core tank (No. 1 in the figure) and the air 


9 


REPRODUCIBILITY OF THSl 
ORIGINAL PAGE IS POOR 



flow guide baffles (No. 2). The reactor was designed to simulate a two 
region gaseous core reactor in which the central fissioning core is sur- 
rounded by a working fluid gas, so the core tank was surrounded by a void 
of 27.14 cm (10.69 in) thickness. n>e outer tank was also constructed of 
aluminum and had an outside diameter of 183.1 cm (72 in), a weight of 170 
Kg (375 lb) and a wall thickness of 0.635 cm (0.25 in). This outer tank 
is shown as item 3 in Fig. II-3. The outer spherical tank was immersed in 
a much larger tank of D 2 O (heavy water) to provide a minimum DgO thickness 
of 91.4 cm (3 feet). Control rod guide tubes (No. 11) extended from the 
side of the tank to permit the insertion of control rods by means of con- 
trol red actuators (J?o. 9). The reactor was started up by turning on the 
electric heaters (J^o. 28) causing uranium hexafluoride vapor to flow through 
the fuel line (No. 14) into the inner tank, so the reactor went critical 
when enough vaporized uranitnn hexafluoride fuel had entered the cavity. 

Prior to start-up, a 19-group transport calculation had been performed to 
determine the critical mass. The calculated critical mass was 8.4 kilograms 
(18.52 lb) of uranium; the actual critical mass turned out to be 8.434 kilo- 
grams (18,59 lb). 

Three configurations were studied in this series of experiments. The 
difference was in the region between the core tank and the cavity tank. In 
the first experiment this void was filled with air, in the second it con- 
tained hydrogen in the form of polystyrene with polyethylene sheets dispersed 
uniformly between the slabs of polystyrene. The third experiment contained 
a stainless steel liner 0.076 cm (0.03 in) thick on the inside of the outer 
tank forming the cavity. The critical mass increased from 8.43 kilograms 
(18,59 lb) with only aluminum baffles in the void region to 29.38 Kg (64.8 lb) 
with both hydrogen and the stainless steel liner added to the void region. 
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Table II-2 provides a comparison of the Soviet and U.S. uranium hexa- 


fluoride reactor experiments. 


Table H-2 Comparison of TJFg Critical Experiments 


Location 

USSR 

Idaho Falls 

Idaho Falls 

Year Operated 

1964-57 

1967 

1969 

Enri chnient 

90% 

93% 

93% 

Moderator 




in core 

Beryllium 

none 

none 

Reflector 

Graphite 

DjO 

DgO 

Core Geometry 

Cylinder 

Cylinder 

Sphere 

Core Diameter 

1.16 m 

1.22 m 

1.26 m 


(3.81 ft) 

(4 ft) 

(4,13 ft) 

Core Height 

1 .08 m 

1.09 m 

1.26 m 

(3.51 ft) 

(3.58 ft) 

(4.13 ft) 

Volume of UFg 

213 1 . 

1274 1 

1047 1 

in core 

(7.52 ft^) 

(45 ft^) 

(37 ft3) 

Reflector 

50 cm 

91 cm 

91 cm 

thickness 

(20 in) 

(36 in) 

(36 in) 

Critical Mass 

2.024 Kg 

21.3 Kg 

8.4 Kg 

of U-235 

(4.46 lb) 

(47 lb) 

(18.5 lb) 


Thermodynamic Cycle Analysis 

The themodynamic cycle analysis in the previous report (Progress Re- 
4 

port No. 1) showed that high efficiencies can be achieved using UFg as the 
working fluid for Ranlcine or Brayton cycles without requiring excessive tem- 
peratures. In both cycles gaseous phase UFg was employed in the reactor 
core, and a regenerative heat exchanger was recommended between the turbine 
and condenser to improve the cycle efficiencies. 
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The plant schematic and T-S diagram for the regenerative Rankine cycle 
are shown in Fig. II--4. The cycle calculations were performed for reactor 
outlet temperatures ranging from 900“K to 1500®K (1160“F to 2240“t) and for 
reactor outlet pressures ranging from 10 to 160 bars (145.04 to 2320.64 psia) . 
The turbine outlet pressure was taken to be 14.6 bars (211^76 psia), which 


corresponds to the condensation temperature at 338 °K (148°F) . The boiler 
feed pump and turbine effeciencies were taken to be 0.88, the condenser pres- 
sure 1,5 bars (21,76 psi) , and the minimum temperature difference in the 
regenerative heat exchanger AT^ was taken to be 28°K (SO^F) . As expected 
the higher pressure and temperatures give the best overall efficiences. For 
a turbine inlet temperature of 1500“K (2700°H) and pressure of 160 bars 
(2320.6 psia), the overall efficiency is 58%. Table 11 - 3 summarizes the 
results of Rankine cycle calculations. 


Table II-3 Rankine Cycle Parameters for - 21.76 psi and T^ - 606. 9®R 




T/R 

n(%) 

W^(Btu/lb ) 
T m 

w .. 

net 

t„”r 

2 


T “R 
3 

P, = 145.04 

psi 








1 


1620 

29.58 

16.54 

16.42 

1458 

1088 

658 



1980 

33.45 

19.88 

19.76 

1788 

1525 

65P 



2340 

37.57 

23.64 

23.53 

2110 

1746 

653 



2700 

41.76 

28.16 

28;04 

2437 

2072 

658 

P, » 580.16 

psi 








X 


1620 

38.02 

24.68 

24.17 

1362 

1027 

66C 



1980 

43.66 

31.05 

30.54 

1668 

132; 

660 



2340 

47.55 

36.23 

35 . 69 

1985 

1634 

660 



2700 

51.98 

43.16 

42.64 

2291 

1937 

660 


= 1450.04 psi 


1620 

41.87 

29.84 

28.52 

1291 

998 

667 

1980 

48.01 

38.07 

36.75 

1591 

1283 

667 

2340 

52.03 

44.49 

43.18 

1898 

1377 

667 

2700 

55.74 

51.18 

49.87 

2199 

1809 

667 


Pj = 2320.64 psi 


1620 

43.40 

32.64 

30,52 

1246 

963 

671 

1980 

50.07 

42.03 

39.92 

1546 

1281 

671 

2340 

54.23 

49,27 

47.16 

1824 

1549 

671 

2700 

58,25 

57,64 

55,53 

2147 

1829 

671 
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ENTROPY ,S 

Fig. II-4 Schematic of Ranklne Cycle Power Plant with Regenerative anc 
its Temperature-Entropy Diagram 
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A plant schematic and temperature entropy diagram for the regenerative 
gas turbine cycle using as the working fluid is illustrated by Fig. II-5. 
Calculations were done for reactor outlet temperatures ranging from 1100“K 
to 1500°K.. Since there is an optimum compressor pressure ratio for the 
maximum efficiency of the Brayton cycle, the turbine inlet pressure was. 
varied from 10 bars to 160 bars (145.04 psia to 2320.6 psia) for fixed com- 
pressor inlet pressure of 1 bar (14.50 psia) and 4.5 bar (65.63 psia). Al- 
so a pressure loss of 0.25 bar (3.63 psia) was assumed for each pass of heat 
exchangers and reactors, because the Brayton cycle efficiency is significant- 
ly depended on the pressure losses in the system. The minimum temperature 
difference between the two fluid streams in the regenerator was set at 28°K 
(50°R) . Turbine and compressor efficiencies were 0.88. Figures II-6 and II-7 
show that Brayton cycle efficiencies as a function of the ratio of turbine 
inlet pressure to compressor inlet pressure. As shown in the figure, the 
maximum efficiency occurred around the pressure ratio of 20. For the same 
inlet temperature the efficiencies for plants with the lower compressor in- 
let pressure of 1 bar (14.50 psia) are higher than the plant with a compres- 
sor inlet pressure of 4.5 bar. This is due to compressor inlet temperatures 
being lower for the lower compressor inlet pressure. 

Comparison of Ranlcine and Brayton cycles shows that for the same reactor 
exit temperature, Rankine cycle can achieve higher thermodynamic efficiency. 
This due to the fact that Rankine cycle reject most of the heat at a low 
temperature than the Brayton cycle. However, this does not mean that Brayton 
cycles are undesirable because their higher temperature for heat rejection 
may make them more suitable for space applications. 

UFg Breeder Reactor Studies 

If the UFg- fueled reactor is to serve as a large-scale energy resource, 
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ENTROPY, S 

Figure II-5 Regenerative Brayton Cycle Schematic and T“S Diagram 
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PRESSURE RATIO P^/p4 

Figure II-7 Regenerative Bray ton, Cycle Efficiency vs. Pressure Ratio 
with P4 = 4.5 bar (65.63 psi) 

either on the ground or in synchronous orbit, it must breed its own fission- 
able fuel. A UF.-fueled fast breeder reactor study was carried out in the 
o 

Soviet Union and reported in 1970.^ The UFg flowed through a large number 
of small tubes joined by manifolds above and below the reactor core, as 
shown in Fig. H-8. Sodium coolant flowing between the tubes removes heat 
from the core. The reactor is brought to criticality at a pressure of 400 
atmospheres. Inconel-X could be used as the tube material since it is cor- 
rosion resistant to UFg and sodium at temperature of around 800°K (IGOO^F) . 
The problem of the radioactive decomposition of the UFg can be eliminated by 



1 

V. 
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Figure II-8 Diagram of a DF^-Fueled Fast Breeder Reactor® 


the addition of a strong fluorinating agent such as Table II-4 provides 

5 

a comparison between this UFg gas fueled fast breeder and the Fermi reactor . 

8 

Table II>-4 Comparison of UFg Fast Breeder with the Fermi Fast Reactor 



Fermi Reactor 

UFg Breeder 

TheriTial Power 

300 MWt 

300 MVJt 

U-235 Critical 
Mass 

485 Kg 
{1069 lbs) 

250 Kg 
(551 lbs) 

Core Diameter 

0.77 m 
(2.53 ft) 

1.01 m 
(3.31 ft) 

Core Height 

0.77 m 
(2.53 ft) 

1.01 m 
(3.31 ft) 

Sodium Inlet 
Temperature 

561 *K 
(550 °F) 

561 '‘K 
(550 °F) 

Sodium (utlet 

700 

(800 °F) 

700 “K 
(800 °Fl 

Velocity of Na 
in Core 

9.9 m/sec 
(32.5 ft/sec) 

7.4 m/sec 
(24.3 ft/sec) 


In the U.S.A. until recently Che nuclear analytical and system studies 
had dealt only with non-breeder UFg reactors. Research was undertaken at 
Georgia Tech to determine the potential for breeding of UFg fueled reactors. 
Tlie 26 energy group, Mach computer code with thermal group cross section, 
prepared by THERMOS, was used to calculate critical masses and breeding ra- 
tios of spherical, UFg gas-core reactors with 70 centimeters (27.6 inches) 
of blankets. The blankets were composed of thorium and beryllium. Ute atom- 
percent of thorium in the blanket was varied from 0.5 to 2.5 and the core di- 
ameter varied from 1 to 5 meters. The average core temperature was taken to 
be 1000 °K. The results of survey calculations showing Che dependence of 
breeding ratio, critical mass, and core pressure upon the atom % of thorium 
in the blanket and core diameter are listed in Table II-5. 
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Table II-5 Breeding Ratio at 1000"K 


Atom % Thorium in Blanket 
Core ■ . 


Uia 

0.5 

1.0 

1.5 

2.0 

2.5 

IM 

1.1229 

1.2161 

1.2581 

1.2913 

1.3093 

2M 

1.1569 

1.2157 

1.2574 

1.2840 

1.3020 

3M 

1.1662 

1.2160 

1.2549 

1.2821 

1.2992 

4M 

1.1755 

1.2166 

1.2536 

1.2795 

1.2977 

5M 

1.1525 

1.2073 

1.2502 

1.2775 

1.2962 


Kg. Uranium 233 Required at 1000 “K 
Atom % Thorium in Blanket 

Core 


Dia 

0.5 

1.0 

1.5 

2.0 

2.5 

ir.i 

30.432 

119.58 

185.15 

223.58 

252.32 


65.612 

368.68 

613.8 

781.31 

898.74 

3m 

122.60 

751.96 

1304,0 

1674.2 

1934.8 


192.20 

1268.9 

2234.4 

2903.3 

2375.4 

5m 

428.95 

2332.9 

3745.8 

4705.9 

5332.2 


iore 

Dia 

0.5 

Core Pressure at 1000 ®K (bar) 

• Atom % Thorium in Blanket 

1.0 1.5 2.0 

2.5 

Im 

20.41 

79.15 

120.0 

150.6 

1?0.1 

2m 

5.90 

30.77 

51.19 

67.9 

75.91 


4,46 

19.24 

32.55 

41.28 

47.78 

4m 

3.51 

14.90 

24,19 

30.15 

35.70 

5m 

3.82 

13.96 

19.99 

26.04 

29.01 
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These results are encouraging since they indicate that a UFg breeder 
reactor with a 1 meter (3.28 ft) diameter and a 120 Kg (265 lbs.) critical 
mass could operate at a core pressure of about 80 bar (1160 psi) with a 
breeding ratio of 1,22, with only 1% thorium in the beryllium blanket. Also 
a reactor with a 2 meter (6.56 feet) diameter and 614 Kg (1354 lbs.) criti- 
cal mass could operate at 50 bar (725 psia) with a breeding ratio of 1.26. 

An analysis of the effect of structural materials on the performance 
of the gaseous core breeder estimated only a few percent loss in breeding 
ratio and less than one percent increase in critical mass and pressure . ^ 

Recent studies on structural materials of the UF^ breeder reactors 
pursued at Georgia Tech indicate that careful selection of coolant in the 
blanket is vitally important in the design of breeder reactor. 

Nuclear Analysis 

Analytical investigations of the nuclear characteristics of gaseous 
core reactors were begun in the mid 1950 *s as interest began to develop in 
the possibility of using reactors of this type for rocket propulsion. These 
calculations were one dimensional in nature and did not take into account 
complicating factors such as temperature distributions within the cavity and 
the effect of up-scattering of neutrons. Figure II-9 illustrates the critical 
mass as a function of reactor diameter and cavity void fraction for a cavity 
with length equal to its diameter, as calculated at the Jet Propulsion Labo- 
ratory using one-dimensional diffusion theory and the Fermi age approximation.^ 
This figure illustrates some of the basic neutronics characteristics of cavity 

reactors. For example, there is a minimum cavity size at which the moderator 

% 

can support a self-sustaining reaction in a thermal reactor. 

At the Douglas Aircraft Corporation calculations were performed for a 
cavity reactor containing a central region of uranium surrounded by an annular 
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Figure II-9 


Pu-239 Critical Mass in a Cylindrical 
BeO Moderated Cavity Reactor^ 
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hydrogen zone contained within the reflector-moderator with a nozzle at the 

end of the chatnher. One and two dimensional transport calculations were 

8 

performed using 14 energy groups . The difference in reactivity determina- 
tions between diffusion and Sg transport theory was only 0.0005. These 
techniques were used to survey a variety of fuel materials, reflector mate- 
rials, core geometries, and reflector thicknesses. Figure II-IO shows the 
critical mass in kilograms of U-235 versus core size for graphite and beryl- 
lium oxide reflectors with data given for various reflector thicknesses. 

For the graphite moderated reactor the minimum critical mass occurred for a 
cavity diameter of about 1.6 m. (5.25 ft.) and for beryllium oxide the mini- 
mum critical mass occurred for a 0.8 m. (2.62 ft) diameter. Beryllium 
oxide is a superior moderator because it slows down the neutrons closer to 
the core so that they have a greater probability of returning to the core 
to cause further fission, rather than leaking out of the system. Beryllium 
moderators typically have more neutron absorption than graphite, whereas 
graphite has a greater loss through leakage. 

The nozzel of the reactor can have a significant effect on critical 
mass, causing it to increase as much as 50%. Both Figs. II-9 and II- 10 fail 
to account for the presence of a nozzle. 

One important neutronics aspect of cavity reactors is that the fuel 
temperature does not reflect very greatly on the neutron temperature. The 
neutron temperature is dominated by the temperature of the reflector. The 
fuel Doppler effect does have a minor influence, but the only important 
effect of the fuel temperature on the system results from the high pressure 
necessary to contain a critical mass of extremely hot fissioning plasma. 

Researchers at the United Aircraft Corporation calculated the critical 
mass requirements for a vortex stabilized gaseous core reactor using one and 
two dimensional diffusion theory.^ The length and diameter of the cylindri- 
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cal cavity was taken to be equal to 1.83 m. (6 ft.), the cavity liner was 
composed of the tungsten-184 isotope (which has a lower cross section than 
natural tungsten) , and the surrounding moderator region consisted of suc- 
cessive layers of beryllium, beryllium oxide, graphite, and heavy water. 

The calculations made allowances for an annular duct leading from the cavity 
to the exhaust nozzle, a fuel-injection port, voids in the moderator, the 
radial distribution of hydrogen temperature in the propellant region, and 
the radial fuel density distribution. In addition the moderator was con- 
sidered to be surrounded by layers of natural tungsten and iron to simulate 
the external piping in the pressure vessel. 

The critical mass of U-233 fuel in the cavity was calculated to be 
22.7 kg. (50.1 lbs.). One dimensional calculations indicated that this 
critical mass could be reduced about 40% by using a cavity liner of beryl- 
lium tubes covered with niobium carbide coated graphite sleeves, instead of 
the ttjngsten-184. This result again indicates the importance of neutron 
absorbing materials at the moderator-core interface. 

Researchers at the Kurchatov Institute for Atomic Energy in Moscow 
investigated acoustical instabilities in boty cylindrical and spherical 
gaseous fueled reactors . The problem was formulated in the age dif- 

fusion approximation by applying the appropriate boundary conditions to 
develop an analytical solution which describes the development and amplifi- 
cation of fuel density oscillations as a function of reactor parameters. 
Acoustical oscillations and also total fuel mass oscillations. In which 
the fuel density varied simultaneously throughout the core volume, were 
considered for a variety of reactor parameters and system geometries. The 
acoustical oscillations were symmetric about the reactor center in the 
spherical reactor and were asymmetric in the cylindrical reactor, that is 
the oscillations propagated from the center to the wall and back again. It 






26 


was shown that the dynamic effects which occur in a gaseous core reactor 
during fuel density oscillations can be completely described by two terms: 
the dynamic reactivity and the variable component of the neutron flux. The 
dynamic reactivity is the net reactivity introduced by an oscillation over 
a complete cycle and the variable component of the neutron flux is the 
amplitude of the oscillation. The dynamic reactivity was shoT'?n to always 
be negative for beryllium and graphite moderated gaseous core reactors, but 
it could be positive for heavy water moderated reactors. Thus fuel density 
oscillations in beryllium and graphite moderated cavity reactors would tend 
to cause the reactor to shut down. Another result was that acoustical oscfl 
lations introduce a negative reactivity effect, except in the case of very 
slow oscillations along the channel axis of a cylindrical reactor, in which 
case the reactivity arises as a first order effect with respect to the 
oscillation amplitude. 

Previous Studies in Georgia Tech (Nuclear analysis) 

Exploratory calculations were performed at Georgia Tech for several 
very high temperature gas core breeder reactor configurations,^ The com- 
putation method involved the use of the MACH-1 one-dimensional diffusion 
theory code^^ and the THERMOS integral transport theory code^^ for thermal 
cross sections. Results of the parametric survey showed that breeding 
ratios in the range of 1.06 - 1.12 could be obtained with criti cal masses of 
300 to 850 kg. of 11-233 for various material compositions in. a 5 meter 
diameter cavity with a 0.5 meter thick blanket. The effects of changes in 
fissile material in the blanket were examined. The addition of structural 
material in the blanket decreased the breeding ratio, by approximately 2% for 
0.2 atom percent natural molybdenum or 4% for fully enriched molybdenum. 

The effect of changes in cavity temperature were also examined, MACtt-1 was 


used as the primary computational tool in these studies, and THEBMOS was 
used to supply thermal neutron parameters to MACH-1. MACH-1 is a one- 
dimensional diffusion theory code with one thermal group (no upscatter) and 
Thermos is a one-dimensional integral transport theory code in the thermal 
range with complete upscattering. All the reactor configurations were 
spherical and hence amenable to one-dimensional analysis. For the MACH-1 
code the 26 group ABBN cross section set of Bondarenko, et al. was used. 

The thermal group of the ABBN set is for 2200 ra/sec. (0.0253 eV) neutrons 
and hence is not realistic for the high temperatures which can be produced 
in a gaseous core reactor. The THEEMOS code was thus used to determine 
thermal cross sections to be inserted into the MACH-1 computation along with 
the ABBN set. For a given configuration the computational method was as 
follows : 

1. Run MACH-1 with 26 group ABBN to estimate critical concentrations 
preliminary results. 

2. Run THERMOS with 50 groups (up to 2.15 eV) using above concentrations 
and calculate spatial and spectral averaged cross sections. 

3. Run MACH-1 with 22 fast groups from ABBN (2.15 eV) and thermal 
cross sections from THEEMOS run. 

Thus the final results of a computation may be thought of as a 23 group cal- 
culation with one thermal group using a thermal cutoff of 2.15 eV. Steps 
2 and 3 were repeated if the final concentrations varied markedly from the 
estimates; steps 1 and 2 were omitted for very similar configurations. The 
high thermal cutoff value is required because of the possibility of a large 
increase in neutron energy due to upscatter in the hot gaseous core. 

Explicit in all calculations were the assumptions associated with the 
two computer codes. Diffusion theory does not seem to be very restrictive 
based on previous comparisons to transport theory for gaseous fueled cavity 
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reactors. For example. Hall found a difference of only 0.0005 in 
between diffusion and Sg transport theory. THERMOS contains the assumption 
of isotropic scattering but this is felt to be reasonable at the energies 
involved (<2.15 eV) . More restrictive assumptions for the THERMOS runs 
were the slowiug-down source and the U-233 resonances below 2.15 eV. The 
slowing-down source was assumed to be spatially independent; MACH estimates 
showed that the epithermal flux is rather flat in the cavity but decreases 
rapidly in the blanket region. This would imply then that the flat source 
assumption is rather good for the cavity and not good in the blanket. Since 
the temperature is not as high in the blanket and resonance capture is 
important in thorium, the results are not as sensitive to thermal cross 
sections for blanket materials as in the cavity. No Doppler broadening 
capabilities exist with THERMOS so these resonances were included at room 
temperature only. The assumption of no Doppler broadening should not be 
too severe since the resonances are very broad even at room temperature. 

Since only eight of the fifty THERMOS groups are used to span these reson- 
ances, the results were probably less sensitive to Doppler broadening than 
to the low number of groups in that interval. 

A general assessment of the computational method and its assumptions 
was provided by analyzing a configuration previously analyzed by Whitraarsh^^ 
with a 3.05 m (10 foot) cavity diameter and 0.61 m (2 foot) thick reflector. 

An essentially equivalent configuration was obtained by reducing the number 
of regions and by homogenizing similar regions, then the computational method 
outlined previously was used to analyze this configuration. The THERMOS 
computational was performed for the cavity regions only. Final results 
gave a value of k of 0,986 as compared with 1.000 for Whitmarsh, a difference 
of only 1.4%. In light of the homogenization used to obtain a nearly equi- 
valent configuration, the agreement indicates that the computational techniques 
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are valid for survey calculations. The largest source of discrepancy was 
attributed to the sensitivity to U-235 thermal cross sections. The agree- 
ment does show that his computational method should be sufficient to identify 
trends and perform parametric studies for various gaseous core nuclear 
reactors. All the configurations examined were spherical in geometry. The 
cavity region contained hydrogen as moderator/coolant and U-233 as fuel. 

The blanket consisted of graphite and thorium. The relative concentrations 
of the materials as well as the size of the regions were varied parametrically 
to examine a matrix of cases in an attempt to obtain the most feasible 
gaseous core thermal breeder concept. Some of the results are illustrated 
in Table II-6. Heat transfer and system analysis studies had estimated the 
bulk average cavity temperature in the reactor core to be in the range of 
3000°K to 4000“K; thermal cross sections and pressures were calculated for 
the case of 4000“K. For the doubling time calculations a power level of 
1000 MWt(th) was assumed. A higher power level shortens the doubling time 
proportionally . 
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Table II-6 Parametric Study of Relative Material Concentrations 
in a Gaseous Core Breeder Reactor^ 


Cavity 
H/U Ratio 

Blanket 
C/Th Ratio 

U-233 Mass 
(Kg) 

Breeding 

Ratio 

H Press 
(atm) 

Doubling 
Tfme (yr) 

140/1 

2/1 

452 

1.1026 

(a) (b) 
710 514 

9.6 


4/1 

390 

1,0962 

612 443 

8.9 


10/1 

301 

1.0636 

472 342 

10.3 

100/1 

2/1 

576 

1.1056 

646 468 

11.9 


4/1 

494 

1.0997 

553 401 

10.8 


10/1 

375 

1.0662 

420 304 

12,4 

.60/1 

2/1 

847 

1.1029 

569 413 

17.8 


4/1 

721 

1.0966 

485 351 

16.2 


10/1 

537 

1.0635 

361 261 

18.4 


(Cavity radius 250 cm, Blanket thickness 50 cm) 

(a) Hydrogen partial pressure at AOOO *K, H 2 mole fraction 0.92, 

(b) Hydrogen partial pressure at 3000 "K, H 2 mole fraction 0,99, 

(c) For 1000 (t), proportionally lower per higher average power. 
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MATERIALS CONSIDERATIONS FOR THE GAS CORE BREEDER REACOOR 

The design of the gas core breeder reactor required selection of 
materials for the reactor fuel, reactor blanket, core liner, reactor pres- 
sure vessel, system piping, and turbine blades. This chapter discusses the 
selection of these materials and summarizes some of the properties needed 
for the system design, 

3.1 Reactor Fuel 

The reactor design criteria specified a gaseous fuel. Since uranium 
hexafluoride is the only chemical compound of uranixim that exists in a gase- 
ous state at reasonable temperatures, UF^ is the selected reactor fuel. 
Appendices A and B list thermodynamic and thermo-physical property data for 

UFg. 

3.2 Reactor Blanket 

Several concepts were considered for the reactor blanket as discussed 
in Chapter VII. The final selection was to use a molten salt for the reactor 
blanket which is similar to the fuel used in the Molten Salt Breeder Reactor 
(MSBR) . 

The fuel salt used in the MSBR is LiF~BeF 2 ~lhF^~UF^ (71.7-16-12-0.3 
mole %). In order to .reduce parasitic neutron capture in the lithium, the 
lithium is enriched to 99.995% in Li^. For use as a blanket for a gas core 
breeder reactor, the MSBR fuel has the UF^ replaced with ThF^ and the blanket 
composition becomes LiF-BeF 2 “ThF^ (71.7-16-12.3 mole %). This salt melts at 


about 930°F, has a low vapor pressure at the operating temperature, and is 
stable in the proposed range of application (IOOO-1600'’F) . Ihe blanket salt, 
with a vsicosity twice that of kerosene, a volumetric heat capacity close to 
that of water, and a thermal conductivity about twice that of water has ade- 
quate heat transfer characteristics and acceptable pressure losses due to 
flow. Some of the properties of the molten salt are given in Appendix C. 

The LiF-BeF 2 -ThF^ salt has been found to be compatible witli the materi- 
als in the system, modified Hastelloy-N, XIFg* etc. In reaching die final 
selection for an adequate blanket salt for the gas core breeder reactor, a 
blanket salt consisting of elements having the following qualities had to 
be picked : 

a) low neutron capture cross section 

b) thermally stable mixture 

c) low operating pressure 

d) non-aggressive toward other materials 

e) able to survive fissioning of U and Pu 

f) Tolerant of fission product acciunulation without serious deteriora- 
tion of its essential properties 

g) low fuel cycle cost dependent on (1) ability to be reprocessed for 
turn around of unbumed fissile material, (2) ability to recover 
bred fissile material, and (3) ability to remove fission product 

poisons . 

1 

Grimes discusses why the fluorides form the only salts with acceptable 
absorption cross sections, stability, and melting temperature. This is due 
primarily to the stability of the ThF^ in the salt since the ThF^ is the 
weak link. Fortunately its high melting temperature is reduced with the 
addition of diluent fluorides (BeF 2 and LiF). 
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For succsssful operation of the reactor, the oxide contamination of the 
blanket must be kept at low levels. In addition, the ratio of LiF to BeF 2 
should be high in order to keep the viscosity at the desirable low level. 

The BeF 2 concentration is advised to be kept in the range of 16 to 25% (for 
a ThF^ concentration of 12%) in order to maintain the liquidus temperature 
below 932 “F. 

3.3 Structural Materials 

3.3.1 Modified Hastelloy-N for Pressure Vessel and Piping 
A material that has been developed for containing the molten salt 
blanket is an alloy which is a slight modification of the present commercial 
Hastelloy-N. The modified Hastelloy-N is very similar in composition and 
other related physical properties to standard Hastelloy-N. However, the 
modified version is superior because of its ability to resist helium embrit- 
tlement due to neutron irradiation. Helium embrittlement is the result of 
helium production in alloys due to tlie (n,Cf) reaction of boron under neutron 
irradiation. Boron is found in minute quantities in most alloys as a result 
of the manufacturing process. However, the amount of boron present is suf- 
ficient to result in serious damage to the alloy when it undergoes neutron 
irradiation . 

The problem of irradiation embrittlement of Hastelloy-N has been solved 

2 

by the addition of 2% titanium to the alloy. Experiments have been per- 
formed on the modified Hastelloy-N in which it has been irradiated with a 
neutron fluence of 10 nvt at 650“C and then stressed at 650“C, It appears 
diat for stresses below 10,000 psi, irradiation has no effect on the stress- 

rupture properties, i.e. for stresses below 10,000 psi the material will 

2 

have an "infinite" life. Figure III-l illustrates stress -rupture charac- 

2 

teristics of modified Hastelloy-N for various neutron fluences. 
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Fig. III-l Stress-Rupture Properties of Several Heats of 

Hastelloy-N Irradiated at 650“G to the Indicated 
Thermal Neutron Fluence and Tested at 650°C. 

(The number by each point indicates the strain.) 
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Hastelloy-N is an alloy developed specifically for use with molten 

fluoride systems* Among the constituents, Cr is the least resistant to 

attack by fluorides. Fortunately, the Cr content of Hastelloy-N is low 

enough for the alloy to have excellent corrosion resistance toward the 

salts. Also, the Cr content is high enough for high oxidizing resistance. 

Corrosion by fluoride salts is by selective chromium removal, thus 

alloys such as stainless steels witli a high Cr content are more susceptible 
2 

to corrosion. For the molten salt used for the gas core breeder reactor 

blanket, the envir^onment is moderately oxidizing and the corrosion process 

is dependent on the diffusion of chromium to the surface of the alloy and 

its subsequent removal. Under these conditions the Hastelloy-N exposed to 

the blanket salt should suffer no greater depth of chromium removal than 

3 

8 to 10 mils over a 30 year Life. If UFg is also only moderately oxidizing, 
then the rate of corrosion will also be dependent on the diffusion of chro- 
mium and hence the depth of chromium removal will also be 8 to 10 mils over 
a 30 year life. 

From the standpoint of the MSBR project, all of the major questions 

with regard to the use of modified Hastelloy-N had been solved with the ex- 

2 

caption of the intergranular cracking of Hastelloy-N by tellurium. The 

problem of intergranular cracking was not "discovered" until late in tlie 

MSBR project when the MSRE (Molten Salt Reactor Experiment) was shut down 

after three years of operation and its components examined. Intergranular 

cracking is the result of the diffusion of tdie fission product tellurium 

along the grain boundaries of the alloy where it combines witli nickel to 

form a weak structural compound. Consequently, when the alloy is then 

2 

stressed, cracks form along the grain boundaries. 
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Tellurium forms a weak fluorine compound which is easily broken down 
resulting in tellurium being deposited on the alloy*s surface thus beginning 
the intergranular cracking process. The cracking process increases with 
increasing tellurium concentration and the rate of deposition of tellurium, 
Ihere is a slight increase in cracking with an increase in temperature from 
650‘’C to 700°C. It has been found that intergranular cracking was signifi* 
cantly reduced and in some cases almost non-existent in alloys with a high 
chromium content, i.e,, 15 to 16% and above. For example, type 3Q4L stain- 
less steel shows a high resistance to tellurium cracking. It is also theo- 
rized that the diffusion rate influences the depth of cracking. This tends 

to explain the effect of chromium, since the chromium concentration is di- 

2 

rectly related to the diffusion rate. Based on conclusions drawn from the 
MSRE analysis, irradiation does not seem to be a major factor in intergranu- 
lar cracking.^ 

In experiments with Hastelloy-N and other alloys, the chromium content 

was seen to have the greatest effect in reducing intergranular cracking, 

although the titanium content also may be beneficial in resisting tellurium 

cracking. This conclusion was drawn mainly from the fact that modified 

llastelloy-N (2% Ti) showed better resistance to cracking than did standard 

Hastelloy-N. With regard to standard Hastelloy-N, minor changes in the 

composition yielded no significant results. Experiments also indicated that 

Al, Nb , and Ce might also help resist cracking. There seems to be some 

confusion in the references as to the effects of stress on intergranular 
2 4 

cracking, ’ A good summary of the results of the study of MSRE components 
(effects of Te) can be found in Ref. 4* 

Since chromium reacts readily with tellurium to form a relatively stable 
compound which doesn*t diffuse readily, it may be desirable to increa^se die 
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chromium content of the modified Hastelloy^N in order to reduce the effects 
of intergranular cracking. However, the question does arise as to whetl^er 
the increased chromium content may cause increased corrosion problems with 
the molten blanket salt. 

In case the intergranular cracking problem of Hastellpy^H cannot be 
solved, type 304L stainless steel is a possible alternative material. This 
is based on the operating experience of a fluoride salt loop built of this 
material. As of August, 1972, this loop had operated satisfactorily for 
9 years at 1250°F and suffered only a 1.5 mil per year corrosion rate (esti- 
mated). However, a close control of the oxidation state is required.. The 
problem of the irradiation embrittlement of the stainless steel has been 
solved by the addition of 0.2% titanium. The remaining question is whether 

or not the corrosion rate can be reduced to an acceptable level by either 

3 

lowering the temperature or by alloy changes. The question then arises 
with regard to alloy changes as to whether or not type 304L stainless steal 
can be made corrosion resistant without reducing its resistance to Te inter- 
granular cracking. 

Another possible alternative is the use of modified Hastelloy-N in 
combination with a Monel 404 liner on the core side. In a telephone conver- 
sation with Jack Davant of Oak Ridge National Laboratory, it was determined 
that Monel 404 could be welded to Hastelloy-N, and both would be corrosion 
resistant. However, at the operating temperatures there could be an expan- 
sion deviation between the two of about 0.25 inch. This number came from 
an extrapolation of the known data at room temperature, since there are no 
published data about Monel 404 at the desired temperatures. 

Theoretically then, modified Hastelloy-N and Monel 404 can be used 
together as reactor materials. Hovrever, if the expansion deviation is much 
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greater than 0.25 inch, this option is removed 


Properties of Monel 404 (Ni-557»; Cu-44%; Mi-0.01%; Fe-0.05%; C-0 
77 “F are given in Table 7X1-1, 

Table III-l 

5 

Properties of Monel 404 at 77“F 

Form and Condition 

Rod, Hot Rolled 

Yield S. (0,2% offset) 

3.1 X 10^ psi 

Tensile S. 

6.9 X lo'^ psi 

Elongation in 2 inches 

45% 

Hardness 


Specific Gravity 

8.88 

Melting Point 

2460 °F 

Tliermal Expansion Coefficient 

7.4 X lO"^ iu./in./“F 

Tensile Mod. of Elasticity 

25 X 10^ psi 


3.3.2 Turbine Blades 

From ANL-7241, it was noted that Inconel alloys had reasonable corrosion 
resistance to UFg as well as good structural properties at 1200°F. Therefore, 
nickel based, high temperature, high strength alloys were checked. Such al- 
loys as T0-Nickel, Alloy 713 C, Alloy 701, Mar M246, M-252, etc. were elimi- 
nated because they either contained boron or cobalt. Boron is undesirable 
because it is a poison which would have to be removed in a reprocessing loop 
before reentering the reactor. Cobalt is easily activated, making turbine 


maintenance difficult. 


Pratt and \^hittvey Aircraft has developed a process for casting high 
strength gas turbine blades from standard nickel based alloys. Tlierefore, 
regular Inconel alloys could be investigated for use in the XIFg turbine. By 
eliminating grain boundaries that lie perpendicular to the principal stress 
axis (this axis lies parallel to the length of the blade), the earlier pro- 
cesses imparted ductility and thermal shock resistance to creep-resistant 
nickel -base super-alloys. The newer process, by eliminating all grain 
boundaries, produces a superior material ^irtien it is applied to Mar-Ml200 
and similar alloys. This is the first of a new family of gas-turbine ma- 
terials that are called "monocrystaloys . " It has, for example, exceptional 
thermal shock resistance since it has not cracked during 2400 shock cycles 

g 

between room temperature and 2200"F, 

The two best Inconel candidates, Inconel 718 and Inconel X-750, were 
examined because of their superior strength properties with respect to all 
otlier Inconel alloys. Based upon a telephone conversation with Luke Yrkovich 
of International Nickel, Inconel X-750 was ultimately chosen because of its 
excellent corrosion resistance to UFg, as well as for its superior strength 
properties . 

Inconel X-750 is an age-hardenable, nickel -chromium alloy used for its 
corrosion and oxidation resistance, and high creep rupture strength at tem- 
peratures up to 1500 “F. The alloy, also containing aluminum and titanium, 
is made age hardenable by the combination of A1 and Ti x»ith Ni to form 
gamma prime — the intermetallic compound Ni^CAljTi). Although much of the 
effect of age hardening is lost with increasing temperatures over 1300°F, 
the fully heat-treated or directly aged alloy will still have useful strength 
up to 1800"F.^ Therefore, it is reasonable to assume that the turbine blades 
should be constructed entirely of Inconel X-750. 


3 .4 Conclusions and Recommendations 


3.4.1 Blanket Salts 

The blanket salt composition was chosen on the basis of neutron cross 
section, viscosity, chemical stability, and liquidus temperature. Ihere is 
little uncertainty with regard to its phase behavior, its physical proper- 
ties, and its interaction with the container and moderator materials. 

g 

More uncertain, however, are the effects of oxidation-reduction state 
of the salt on its surface tension, and on the behavior of the noble metal 
fission products. Significant limitations on the blanket salt are imposed 
by the rather high liquidus temperature (932 °F), and the restricted choice 
of containment material. 

The problem that should be looked into is the production of relatively 
large amounts of tritium by neutron interaction with the lithium. 

Variations of the blanket composition are possible and may prove desir- 
able, to mitigate some of the aforementioned limitations. The concen- 

tration may be varied as required, even for such things as optimization of 
the breeding ratio. The processing is expected to keep the oxide concentra- 
tion in the fuel low. 

The constraints of high LiF temperature and high tritium levels cannot 
be circumvented. Since the reactor is to breed with thermal neutrons, the 
cross sections limit the choice of diluent salt constituents to the fluo- 
rides of Be .and Li. The tritium production could be cut to not much more 
than fission yield, if a NaF-ZrF^-UiF^ blanket salt were used. But, neutron 
absorption in the Na and Zr would eliminate breeding. If the reactor is to 
breed, there is no alternative. 

3.4.2 Component Ma terial s 

Hastelloy-N, modified with 2% Ti, is the main choice for the core and 
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blanket construction. Although additional work is required to solve the 
problem of telluritim intergranular cracking, preliminary results indicate 
this problem can be resolved satisfactorily. 

With regard to the design stresses of modified Hastelloy-N, if extrapo- 
lation of values from the test results does not prove to be satisfactory, 
then the operating temperature and total neutron fluence that the components 
are subject to must be reduced. Theoretically, the irradiation should be 
limited to the extent that the creep ductility will not be less than 5 %. 
Based on the available data, a maximum design stress of 10,000 psi is 

recommended when the material is subject to 650 “C and a neutron flucnce of 
20 

10 nvt over a 30 year period. 

Inconel X-750 is recommended for use as the turbine blade material in 
the po?7er conversion unit. 
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One and Two Dimensional Diffusion Theory Calculations 
12 3 

Several analyses ’ ’ of the gas core nuclear reactor have been per- 
formed, For examining a broad range of designs, one may utilize multigroup 
diffusion theory to observe trends and perform parametric studies in order 
to identify concepts for further study. 

The fii-i.t phase in performing exploratory nuclear analysis for the gaseous 
core nuclear reactor involved implememting the necessary computational tools 
and formalizing a computational method. The major portion of the effort 
early in the study was devoted to this area. In order to expedite this 

4 

phase, the MACH-I code was used as the primary computational tool in the 

nuclear analysis. To allow a more realistic model of theirmal neutron 

processes in the high temperature gaseous core reactor concept, the THERMOS 
5 

code was implemented to supply thermal neutron parameters to MACH-I. 

One dimensional survey calculations were carried out for a UF fueled 
core surrounded by a molten salt blanket as shown in Fig. IV-1. Table IV-1 
lists the composition of the coolant in the blanket. Later, two-dimensional 
calculations were made with r-z geometry using the EXTERMINATOR program. 

The spherical reactor configurations were amenable to one-dimensional 
analysis. For the MACH-I code the 26 group ABBN cross section set of Bon- 
darenko, et al.^ was used. The thermal group of the ABBN set is for 
2200 m/sec (0.0253 eV) neutrons and hence is not realistic for the high 
temperatures which can be produced in a gaseous core reactor. In some cases 
the thermal cross sections were taken as Maxwellian. In other cases the 
THERMOS code was used. The THERMOS code was thus used to determine thermal 
cross sections to be inserted into the MACH— I computation along with the 
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TABLE IV~1 


COMPOSITION OF COOLANT IN BLANKET 



MOLE % 

VOLUME % 

DENSITY (Gr/cm^) 

* 




LiF 

71.7 

42.28 

2.635 

BeF^ 

16.0 

22.08 

1.986 

ThF, 

12.3 

35.00 

6.32 

* 

7 

6 



7 6 

The Li is 99.99% Li because Li has a large thermal absorption 


cross section. 




ABBN set. For a given configuration the computational method was as follows: 

1. Run MACH-I with 26 group ANNB to estimate critical concentrations 
preliminary results . 

2. Run THEBHOS with 50 groups (up to 2.15 eV) using above concen- 
trations and calculate spatial and spectral averaged cross sections. 

3. Run MAGH-I with 22 fast groups from AKNB (2.15 eV) and thermal 
cross sections from THERMOS run. 

Thus the final results of a computation may be thought of as a 23 group 
calculation with one thermal group using a thermal cutoff of 2.15 eV. 

Table IV-2 lists the energy groups for the MACH Calculations. 

Steps 2 and 3 were repeated if the final concentrations varied markedly 
from the estimates; steps 1 and 2 were omitted for very similar configura- 
tions. The high thermal cutoff value is required because of the possibi- 
lity of a large increase in neutron energy due to upscatter in the hot 
gaseous core. 

Implicit in all calculations were the assumptions associated with 
the two computer codes. Diffusion theory does not seem to be very re- 
strictive based on previous comparisons to transport theory for gaseous 
fueled cavity reactors. For example. Hall found a difference of only 
0.0005 in between diffusion and Sg transport theory. THERMOS contains 

the assumption of isotropic scattering but this is felt to be reasonable 
at the energies involved (<2.15 eV) . More restrictive assumptions for the 
THERMOS runs were the slowing-doim source and the U-233 resonances below 
2.15 eV. The slowing-down source was assumed to be spatially independent; 
MACH estimates showed that the epithermal flux is rather flat in the cavity 
but decreases rapidly in the blanket region. This would imply then that the 
flat source assumption is rather good for the cavity and not good in the 
blanket. Since the temperature is not as high in the blanket and resonance 
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TABLE IV-2 


MACH-I ENERGY GROUP SCHEME 


GROUP NO. 

ENERGY 

RANGE 

1 

6.5-10.5 

MeV 

2 

4. 0-6. 5 

MeV 

3 

2. 5-4.0 

MeV 

4 

1.4-2. 5 

MeV 

5 

0.8-14 

MeV 

6 

0.4-0. 8 

MeV 

7 

0.2-0. 4 

MeV 

8 

0.1-0. 2 

MeV 

9 

46.5-100 

keV 

10 

21.5-46.5 

keV 

11 

10.0-21.5 

keV 

12 

4.65-10.0 

keV 

13 

2.15-4.65 

keV 

14 

1.0-2.15 

keV 

15 

465-1000 

eV 

16 

215-465 

eV 

17 

100-215 

eV 

18 

46.5-100 

eV 

19 

21.5-46.5 

eV 

20 

10.0-21.5 

eV 

21 

4.65-10.0 

eV 

22 

2.15-4.65 

eV 

23 

1.0-2.15 

eV 

24 

0.465-1.0 

eV 

25 

0.215-0.465 

eV 



capture is important in thorium, the results are not as sensitive to 
thermal cross sections for blanket materials as in the cavity. No Doppler 
broadening capabilities exist with THERMOS so these resonances were included 
at room temperature only. The assumption of no Doppler broadening should not 
be too severe since the resonances are very broad even at room temperature. 
Since only eight of the fifty THERMOS groups are used to span these reso- 
nances, the results were probably less sensitive to Doppler broadening than 
to the low number of groups in that interval. 

A general assessment of the computational method and its assumptions 
was provided by analyzing a configuration previously analyzed by Whitmarsh 
with a 3.05 m (10 foot) cavity diameter and 0.61 m (2 foot) thick reflector. 
An essentially equivalent configuration was obtained by reducing the number 
of regions and by homogenizing similar regions, then the computational 
method outlined previously was used to analyze this configuration. The 
THERMOS computational was performed for the cavity regions only . Final 
results gave a value of k ^^ ^ of 0.986 as compared with 1.000 for Wliit- 
marsh, a difference of only 1.4%. In light of the homogenization used to 
obtain a nearly equivalent configuration, the agreement indicates that the 
Computational techniques are valid for survey calculations. 

The major goal after completing the one-dimensional calculations was 
to do two-dimensional calculations. Instead of using a spherical core as 
was done previously using MACH-I, a cylindrical core was analyzed utilizing 

7 

EXTERMINATOR II, a two-dimensional diffusion theory code. This code 
enabled us to get fluxes in the (r-z) direction and, from these point- 
wise fluxes, the reactor breeding ratio could be calculated. 

The first problem encountered when using EXTERMINATOR II is how to 
calculate the region and energy dependent microscopic cross-sections which 
are input for EXTERMINATOR. This was achieved by using MACH-I for 26 
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groups and collapsing these to 4 groups (thermal) for use in EXTERMINATOR 
The MACH-I cross-section library, however, did not contain any cross 


sections for fluorine and these cross sections (a , a ) had to be found from 

a s 

another source and then input into the MACH-I library before being collapsed 
into 4 groups. The cross-sections obtained for fluorine were from a 123 
group set and had to be collapsed into a 26 group set for use in the MACH 
library. The 26 group flux spectrum from a previous MACH-I run was used as 
the weighting parameter when collapsing the 123 group set into 26 groups. 

Tlie collapsed fluorine cross sections were then input into the MACH-I 
cross-section library where 4 group cross-sections were obtained for the 
desired dimensions and compositions. Dashed lines on Table IV-1 indicate 
the four energy groups for EXTERMINATOR. These region and energy depen- 
dent cross-sections were then used in EXTERMINATOR. Various EXTERtllNATOR 
runs were made to get equal to 1.0, with the dimensions of the core 

being the varied parameter to get equal to 1.0. The height of the 

core was taken to be twice the radius. The core was surrounded by a half- 
inch Hastelloy-N liner and this was in turn surrounded by the molten salt . 
blanket. This was the same arrangement as in the previous one-dimensional 
case. 


Figure IV-2 is a diagram of the geometry used for the two-dimensional 
EXTERMINATOR calculations. Three axial regions were used in the core to 
match the axial temperature distribution more accurately. Table IV-3 sum- 
marizes assumptions made in the calculations. 






TABLE IV-3 


TWO DIMENSIONAL CALCULATIONS 


— CYLINDRICAL CORE WITH REFLECTOR WAS DIVIDED INTO 3 REGIONS 
AXIALLY. 


— TEMPERATURE DISTRIBUTION IN THE CORE IS LINEAR WITH Tin = 
1071°R, Tout = 1560°R 

— AVERAGE TEMPERATURE IN EACH REGION WAS USED TO COMPUTE INPUT 
CROSS SECTIONS 

— 26 ENERGY GROUP INDIVIDUAL ISOTOPE INPUT CROSS SECTIONS WERE 
COLLAPSED INTO REGION-AVERAGED 4 ENERGY GROUP CROSS SECTIONS 
USING MACH-I COMPUTER CODE 


— 4 GROUP CROSS SECTIONS FROM MACH-I WAS INPUTED TO EXTERMINATOR- I I 

TO CALCULATE 

(1) CRITICAL DIMENSION 

(2) CRITICAL MASS 

(3) FISSION DENSITY AS FUNCTION OF POSITION 


IN THE CALCULATION OF BREEDING RATIO, 26 GROUP CROSS SECTIONS OF 


Th 


232 


AND U^^^ WERE COLLAPSED INTO REGION-AVERAGED 4 ENERGY GROUP 


CROSS SECTION MANUALLY WITH FLUX WEIGHT. 
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Heat Generation in Molten Salt Blanket 

Neutrons give rise to heat generation, in non-fissile materials through 
four primary interactions; elastic scattering (n,n); charged particle 
ejections, particularly (n,a) and (n,p) ; inelastic scattering (n,n’y) j 
radiative capture (n,Y) • These four processes may be conveniently divided 
into two classes, direct and indirect, according to whether the heat is 
deposited at the site of the primary interaction or some distance away. 

Elastic scattering and charged particle ejection reactions are direct 
heat generation processes in that the energy of the struck nucleus, or the 
alpha particle or proton, is converted to heat in the solid or liquid within 
distances of less than a millimeter from the site of the primary interaction. 

Radiative capture is an indirect heat generation process since the 
gamma rays that are produced in the reaction are absorbed at various dis- 
tances, often many centimeters, from the site of the primary interaction. 

The primary interaction, therefore, serves only as a source for the gamma 
rays. The kinetic energy given the nucleus as it recoils to conserve momen- 
tum when emitting the gammas is generally only of the order of 10 eV, and so 
may be neglected in heating calculations. 

Inelastic scattering gives rise to both direct and indirect heat 

i 

generation. The daughter nucleus, which recoils to conserve momentum upon 
emitting the inelastic neutron, deposits its kinetic energy at the site 
of the primary interaction. The gamma rays which are subsequently emitted 
by the excited nucleus deposit their energies at various distances frOm the 
site, just as do capture gammas. 

The following energies per core fission were estimated deposited in 
the blanket; 
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TABLE IV-4 


ENERGY DEPOSITION IN BLANKET PER FISSION 


Reaction 


Neutrons from U-233 Fission 
Prompt Gammas 
Fission Product Gammas 
(n,Y) in Thorium-232 
Fissioning in Blanket 
B” from Tli-233 Decay 
3 ” from g^Pa^^3 Decay 


Energy Deposition/Fission 
(MeV/Fission) 

k 

3 

0.5 

9 

1.96 

0.5 

0.1 


Total 19.06 


Because at the present time all calculations have been done for 
spherical geometries it was decided that a rough calculation of heat de- 
posited in the blanket was sufficient. A more detailed analysis should 
wait until a two dimensional code can be used to calculate neutron fluxes 
and reaction rates. 

In the fission process of U-233, 5 MeV goes into kinetic energy of 
the neutrons. Since almost all of the thermalization is done in the blanket, 
it was estimated that 4 MeV was deposited there. It was assumed that 3 
MeV of the 7 MeV prompt gamma energy is deposited in the blanket. There 
will be some gamma energy deposited from fission product decay. This 
is normally 7 MeV; however, since the core is continuously being cleaned 
by reprocessing, most of these decays occur away from the core and it was 
assumed that only .5 MeV is deposited in the blanket. 

When Th-232 absorbs a neutron, it undergoes an (n,y) reaction to reach 
the ground state. The gamma energy is 7.5 MeV and since there are 1.2 reactions 
per fission with the breeding ration currently calculated, 9 MeV is de- 
posited from this source. From the MACH code it was found that 1% of the 
total fissions occurs in the blanket. Thus 1.96 MeV is available per fission 
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in the system. The thorium decay scheme has two decays: 


233p^ 233„ 

These B”’s have energy of 1.23 MeV and .25 MeV. Multiplying by 1.2 gives 
1.5 MeV and .30 MeV per fission in the system. 

Adding these energies gives a total of 20.25 MeV/fission. Since there 
is 196 MeV of recoverable energy per fission, the percent of heat generated 
in the blanket is about 9.72% 

There is a Li^(n,y) reaction, but the atom density of Li^ is so low 
that this energy is negligible compared to that from other sources. 


Results and Conclusions 


Figure IV-3 presents a graph of the one dimensional results showing the 
effect of blanket thickness on the breeding ratio and critical mass. A 
blanket thickness of 90 centimeters corresponds to a breeding i ,io of 1.18, and 
the breeding ratio does not increase as more blanket is added. On concludes 
that 90 centimeters is an optimum choice for the blanket thickness. 

Figures IV-4 and IV-5 presents radial and axial power density distri- 
butions due to fissioning from the EXTERMINATOR results. Notice the fission 
density is relatively flat in the radial direction. The discontinuity in the 
axial fission density (Figure IV-5) is of course due to the regionwise represen- 
tation in the EXTERMINATOR input, and is not a real phenomenon. 

Figure IV-6 is a graph of the energy dependent fission density in the core. 
The median fission energy is 300 keV. Figures IV-7 and IV-8 are the group 
fluxes vs energy averaged over the core and reflector, respectively. The median 
flux in the core is 600 keV in the core and 34 keV in the blanket. 

Figures IV-9 and IV-10 are graphs which show the effect on the spectrum 
of adding moderator (Be and graphite, respectively) to the molten salt blanket. 
As shown in Figure IV-9, the median neutron energy increases from 34 keV (zero 
% Be) to 73 keV (25% Be) and drops to 4.65 eV (75% Be). For carbon the median 
energy increases to 56 keV (25% C) and drops to 7 keV (75% C) . 

When changing the composition of blanket material from 100% molten salt 
(LiF, BeF 2 and ThF^) to 75% molten salt and 25% carbon, we are reducing the 
atom density of lithium in blanket, which is a better moderator than carbon; 
therefore, the spectrum in the blanket region becomes harder. I^hen the volume 
percent of carbon increases to 50%, there is a large increasing of atom density 
of carbon and a small decreasing of atom density oi lithium. Since the ef- 
fectiveness of moderator has to do with both atom density and cross sections, 
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the increasing of atom density of carbon may overcome the fact that some 
lithium is replaced by carbon (or beryllium) . Table IV-5 and IV-6 show the 
effect on median energy, breeding ratio and critical radius when the composi- 
tion of the blanket material is changed from 100% molten salt to 100% carbon 
and beryllium, respectively. 

Table IV-7 summarizes and compares results of the one dimensional 
MACH calculations and the two dimensional EXTEBMINATOR calculations. The 
four group spherical case was in reasonable agreement with the 26 group case, 
so one concludes that the change in going to four groups had a negligible 
effect on the calculated breeding ratio. The breeding ratio of 1.21 for 
the two dimensional reactor looks promising. 

F urther Work 

Many areas, especially those itemized in Table IV-8, must be studied 
further in order to achieve a complete nuclear analysis. 


BLANKET 


THICKNESS (cm) 


FIG. IV-3. MACH-I SURVEY OF BLANKET THICKNESS EFFECTS ON BREEDING RATIOS AND CORE RADIUS 












FLUX INTEGRALS 







Median Energy (keV) 



FIQ. IV-10. GRAPH OF VOLUME % OF CARBON IN BLANKET VS MEDIAN NEUTRON ENERGY 
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PERCENT OF CARBON 
IN BLANKET 

0 % 

25% 

50% 

75% 

100% 

MEDIAN ENERGY (MeV) 
IN CORE 

.53147 

.55617 

.57032 

.57663 

.53016 

MEDIAN ENERGY (keV) 
IN BLANKET 

■34 

46.5 

34 

7.33 

.0252 X 10”^ 

BREEDING 

RATIO 

1.1827 

1.196 

1.190 

1.1331 

0 

CRITICAL 
RADIUS (cm) 

58.5916 

60.9031 

62.6357 

61.3982 

39.2433 ' 

CRITICAL 
HASS (kg U-233) 

379.13 

385.80 

463.22 

436.30 

113.55 







TABLE IV-6 


CRITICAL PARAMETERS VS VOLUME PERCENT OF BERYLLIUM IN BLANKET 


PERCENT OF Be 
IN BLANKET 

0 % 

25% 

50% 

75% 

100% 

MEDIAN ENERGY (MeV) 
IN CORE 

,03147 

.55239 

.55123 

.52800 

44387 

MEDIAN ENERGY (keV) 
IN BLANKET 

34 

73.3 

34 

4.65 X 10"^ . 

0252 X 10"^ 

BREEDING 

RATIO 

1.1827 

1.2284 

1.2032 

1.0651 

0 

CRITICAL 

RADIUS 

58.5916 

61,8469 

61.1468 

53.4241 

29.80 

CRITICAL 
HASS (kg U-233) 

379.13 

445.94 

430.96 

287.43 

49.88 


BLANKET THICKNESS = 114.38 cm 






TABLE IV- 7 

COMPARISION OF CRITICAL PARAMETERS 


o 


BREEDING 
RATIO 


CRITICAL 
RADIUS 


ICRITICAL CORE 
VOLUME 


SPHERICAL CORE 
(26 GROUP) 


1.181 


58.59 cm 


8.425 X loW 


CRITICAL 

HASS (kg U-233) 


379 


BLANKET 
THICKNESS 


114.38 cm 


SPHERICAL CORE 
(4 GROUP) 


1,179 


60.91 cm 


9.466 X 10 V 


426 


114.38 cm 


CYLINDRICAL CORE 
(4 GROUP) 

1.219 

54.78 cm 

1.033 X loW 


496 



Table IV- 8 Neutronic Areas for Further Study 


1. Study more thermal reactor with additional moderation in blanket 

Effect on spectrum 
Effect on breeding ratio 
Effect on control 

2. Select control mechanism 

Reactivity effectiveness of control mechanism 
Effect on breeding ratio 
Effect on power distribution 

3. Calculate reactivity effects 

Temperature and other reactivity coefficients 
Effect of structure and impurities 
Reactivity effects with lifetime 

4. Iterations with thermal and fluid flow analysis 

5. Heating rate in blanket 

6. Dynamic and stability analysis 

7. Analyze more check point experiments - improve model 

8. Potential of co-axial flow gas core reactor - MHD 

9. Investigate applications, especially use of gas core reactor for 

actinide burnup 
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V. 1HERtK)DYWAMIC CYCLE ANALYSIS 

In a previous report^ the use of both Brayton and Rankine cycles 
employing UFg as a working fluid was evaluated. This report showed that 
Rankine cycles gave higher overall plant efficiencies. However, Brayton 
cycles should not be considered completely undesirable because their higher 
temperatures for heat rejection may make them more suitable for space ap- 
plications. However, since this design is concerned with land-based power 
plants, the thermodynamic analysis was conducted using only Rankine cycles. 
The analysis of Ref. 1 is extended to include the evaluation of cycles with 
and without reheat capabilities. In. addition, temperatures and pressures 
at certain selected locations were varied to optimize plant efficiency. 

5.1 Cycle Evaluation 

Calculations were performed for the Rankine cycle illustrated in Fig. 
V-1. Calculations were also done for a cycle with no reheater to determine 
how the cycle efficiency is affected by the reheater. 

Because material limitations of the duct walls limited UFg temperature 
to less than 1660‘’R (921.89°K), the average temperature at the outlet of 
the reactor was chosen to be 1560°R (866.33*’K). The reactor outlet pressure 
was set at 1450 psia (99.97 bars) which is approximately the pressure re- 
quired from the core physics calculations to obtain a critical reactor. 

The minimiim temperature difference in the regenerative heat exchanger, 
4Tg_^, was specified as 50°R (28 “K). The boiler feed pump and turbine ef- 
ficiencies were taken to be 0.88 and the condenser pressure was maintained 
at 21.76 psia (1.5 bars). The pressure loss across the reactor was assumed 
to be 14.89 psi (1,03 bars) and the pressure loss across each side of the 
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regenerator and reheater was taken to be 7.445 psi (.51 bars). In these 
calculations, the pressure loss through the pipes was neglected because the 
plant design was not complete and therefore the length of the pipes between 
components was unknown. 

5.2 Cycle Parameter Calculations 

The thermodynamic cycle calculations were performed using the method 
outlined in Ref. 1. The thermodynamic properties of UF^ are tabulated in 
Appendix A. 

Initially, cycle parameters were calculated without reheating. Ihis 
gave a plant efficiency of 39.03%. Then reheat was added between stages of 
the turbine using the heat generated in the blanket. It was apparent that 
an optimum pressure existed for removing the UFg from the high pressure tur- 
bine and, consequently, P 2 was varied. Figure V-2 illustrates overall plant 
efficiencies as a function of the high pressure turbine outlet pressure. 

The maximum efficiency was 41.44% at a pressure of 435.12 psia (30 bars). 

From Fig, V-2 it appears that the efficiency begins to increase for 
pressures greater than 720 psia (50 bars). However, higher pressures cannot 
be used because the calculations indicated that for these higher pressures 
the reheat temperatures exceed 1560 °R (866 °K), which was not allowed because 
of material limitations. 

Table V-1 summarizes the parameters for the optimized UFg gas-core 
reactor power plant. 



Overall Efficiency = 41.44% 


Tl^ , . =71 = 0.88 

turbine pump 

Pressure Losses = 1 bar in reactor and 1/2 bar in each heat exchanger pass 


State Points tor Stream Locations Shown in Fig, V-1 

#1 B. 


T = 1560. 0“R (866.33“K) 

P = 1450.00 psi (99.97 bars) 
h = 123.88 Btu/lb (288.14 KJ/Kg) 


T = 1516.3 “R (842.06°K) 

P = 427.67 psi (29.49 bars) 
h = 121.63 Btu/lb (282.91 KJ/Kg) 


i5 

T = 1036. 8 “R (575. 65 “K) 

P = 1464.89 psi (100.99 bars) 

h = 62.17 Btu/lb (144.61 KJ/Kg) 
m 


T = 617.3 “R (342. 58 "K) 

P = 1472.34 psi (101.51 bars) 
h = 35.26 Btu/lb (82.01 KJ/Kg) 


#9 


T = 1455.5 “R (803.26“K) 

P = 435.12 psi (30.00 bars) 

h = 114.99 Btu/lb (267.46 KJ/Kg) 
tn 


T = 1293. 8°R (718.46“K) 

P = 29.21 psi (2.01 bars) 

h = 100.86 Btu/lb (234.60 KJ/Kg) 
m 


B. 

T = 667.3 “R (370.36 “K) 

P = 21,76 psi (1.50 bars) 

h = 40.03 Btu/lb (93.09 KJ/Kg) 
m 


£8 

T = 606. 9 °R (336. 83 “K) 

P = 21.76 psi (1.56 bars) 
h = 35.26 Btu/lb (82.02 KJ/Kg) 


T = 606. 9°R (336.83®K) 

P = 21.76 psi (1.56 bars) 
h = 0.00 Btu/lb„ (0 KJ/Kg) 
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VI. HEAT TRANSFER AND FLUID FIX3W 


Because of the high power density in the gas core reactor, it is 
necessary to analyze the heat transfer and flow characteristics of the core 
in order to assure that unacceptably high temperatures are not achieved in 
the core. This requires solving the energy equation for the UFg flowing 
through the core. Because of radial symmetry in the core, the energy 
equation is given by 


pCpU^(r,z) 


^ i. 

bz " r £>r \ 


r K 


E br/ 


+ q'" Cr,z) 


( 1 ) 


where 

u (r,z) = axial velocity 
2 # 

p = density 

= specific heat at constant pressure 

'‘e " + I“=p 'h 

q'" ss volumetric heat j,eneration r-ate 
s„ = eddy diffusivity for heat transfer 
Equation 1 is extremely complex because the UFg physical properties are 
highly temperature dependent and the volumetric heat generation term is 
spatially dependent due to the variable density and neutron flux dis^ 
tributions. The thermal -physical properties of UEg are given in Appendix B. 


6.1 Problem Formulation 

Equation 1 was solved for two sets of boundary conditions: (i) an 

insulated liner wall in which no heat crosses the wall and (ii) insulated 
Liner wall until the wall temperature reaches 1660^11 and then a constant 


wall temperature set at 1660 °R for the rest of the core length. Mathe- 
matically the boundary conditions are given by: 

Case 1 

T(rjO) constant (2a) 


br| 

r=r 

o 


0 for 0 S z 


L 



Case 2 


T(r,0) — constant : (3a) 


with 

= 0 if T g 1600“R 
br w 

(3b) 


<r=r^ ■ 
■0 


and 


(3c) 


T = 1660 °R once T = 1660“R 
w w 


where r ia the core radius and L is the core length. 

■ :■ O ' ■ ■ ■ 

The axial velocity distribution was assumed to be given by a 1/7-power 


laWj or 


u (r,z) — u „„ (1-r/r ) 

z ’ z. Max o' 


1/7 


'■ 0.8 


• O' 


(4) 


•Where u ,, is the maximum (centerline) fluid velocity and u is the 

Z, Max 2, iierui. 

channel average velocity for a given, axial location^ Since the fluid is 
accelerating down the channel length because of the decreasing density 
due to heating, the velocity: distribution in the axial direction was modi- 
fied by multiplying Eq. 4 by the ratio of the average inlet density to the 
a-yerage fl'uid density at a given axial location* Therefore the velocity 
distribution incorporated in E4. 1 is 


( 5 ) 






Ti^here u, is the mean velocity at the core inlet* 

Xj Mean 

In order to calculate heat transfer within a fluid in a forced coa- 
yection (turbulent) system, a property to explain the thermal mixing is 
needed. In turbulent flow the fluid has both radial and axial velocity 
components (the time-averaged radial component is zero) and the mixing of 
two regions ; (caused by turbulent convection) of different temperatures must 
be accounted for. The term to account for this mixing is the eddy diffus- 
ivlty for heat transfer, e^. A basic assumption itl dealing With is that 
it is approximately equal to the eddy diffusivity of momentum transfer, 

(this holds for all fluids except liquid metals) . Numerous expressions 
exist for These expressions are normally obtained from differentiating 

einpirical equations for velocity distributions. One of these expressions 
which is widely accepted is that of Dwyer Which yields a value of zero at 
the channel centerline and wall. The expression for e is in the following 
form^ : 

= (n/p) .72(1/8 £)^ R(L-R^) (6) 

where 

[i = viscosity 
p = density 

£ = Moody friction factor 

= Reynolds number ■ 

E = r/ro 

r = radius to wall 
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The radial dependence of the volumetric heat generation rate was 
modeled from data obtained from a MAGH-I core physics calculation of the 
power density in a spherical gas core reactor. In the axial direction the 
volumetric heat generation, rate was assumed to decrease linearly to one- 
third that at the core entrance. In mathematical terms q'" is given by 

r -0.19(r -r)-jp 2zl 

= q*;*|l + 1.07 e JL^ - 3L J 

Equation 1 was solved numerically by using finite difference representa 
tions for the partial derivatives and incorporating the terms given by Eqs. 

4 to 7. A marching technique was employed which required iteration at each 
axial step in order to incorporate the temperature dependence of the TIFg 
physical properties. The numerical methods used for solving Eq. 1 are given 
in Appendix D. 

6.2 Results of Heat Transfer Calculations 

It was estimated that 9.7% of the reactor power would be deposited in 
the blanket. Consequently, for a reactor power level of 1000 Mwth, 903 Mwth 
would be generated in the reactor core. From the thermodynamic analysis 
described in Chapter V it was found that the core inlet temperature would 
be 1036.8°R (575, 7®K). For the specified outlet temperature of 1S60“F 
(866, 3°K), the UFg mass flow rate through the core must be 13,900 Ib^^^/sec. 

The selected core geometry was a right cylinder with a 100 cm radius 
and a 200 cm length, Flow through the core entrance and exit is through 
numerous inlet and exit nozzles which permits assuming the velocity distri- 
bution given by Eq, 5, Since the core power is 903 Mw, integration of Eq. 7 
over the core volume gives the cons tant qj*' , 


Figure VI-1 illustrates core linear wall temperatures and UFg fuel 
temperatures at the core axis as a function, of core length for the jhsulated 
wall boundary condition given by Eq. 2, As shown, by the calculations, after 
50 cm down the channel length the liner wall temperature exceeds 1660 ?R 
which is considered unacceptably high. Figure VI-2 illustrates the radial 
dependence of UP g temperatures for various axial positions. The tempera- 
tures reach a peak at the wall because of the insulated wall, the volu- 
metric heat generation rate has a maximuia at the wall, and in particular 
because the fluid velocity at the wall is zero ^ich means energy is trahs- 
f erred at that location only ' through conduction. 

Figure VI-3 illustrates core liner wall tenq)eratures and UFg fuel tem- 
peratures at the core axis as a function of core length for the boundary 
condition that the liner wall not exceed 1660 °R (boundary Condition given 
by Eq. 3). Figure VI-4 illustrates the radial dependence of UFg temperatures 
for various axial locations. The UFg temperature at a given axial location 
has a maximum near the channel wall because the fluid velocity is zero at 
the wall. The maximum UFg temperature occurs at the channel exit and is 
2200*’R ^ich is far below temperatures required for substantial llFg ionisa- 
tion. Figure B-1 of Appendix B illustrates UFg composition as a function 
of temperature. 

The boundary condition that the liner wall temperature hot exceed 
1660‘*R requires wall cooling after about 40 cm down the cOre length. Conse- 
quently, it is necessary to examine wall heat fluxes in order to determine 
the extent of the wall cooling. Figure VI- 5 illustrates liner wall heat 
fluxes as a function of channel length. The maximum heat flux occurs at 
the channel exit and has a value of 650 Btu/hr-ft (0.205 watts/cm ). This 
is a small heat flux for which it would be easy to provide wall cooling. 











AXIAL DISTASCE (cm) 


EIS. Vl-3. WMiL XE^IPERATORE VS AXIAL DISTMCE (CASE 2) 











REPRODTJGIBILITZ OF 'THE 
ORIGINAL FAGS IS POOR 








The low thermal conductivity of UFg accounts for the low heat flux in 
the liner wall. The total wall cooling is less than 150 kw which repre 
sents a negligihle amount of the core power. 
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7.1 Design Criteria 

The basis for this project was chosen to be a 1000 MW thermal, station- 
ary power plant employing IIFg-fueled gas-core breeder reactor. The blanket 
for the reactor was initially chosen as a mixture of beryllium and thorium. 

Based on this brief description the actual reactor could be designed 
and modified to take advantage of new ideas and information. The selected 
reactor configurations is shown in Fig. VII-1. This configuration demon- 
strates that a practical UFg gas-core breeder reactor can be constructed 
and this design will serve as a good starting point for future design work 
on this concept. 

7.2 Preliminary Design Concepts 

Based on the design criterion, a multitude of possible designs can be 
described. Numerous combinations were examined and the following concepts 
were considered the most promising or practical to be studied in greater 
depth-: ■ ^ ■ 

-1* A gaseous DFg core with a solid matrix blanket. (Exterrial blanket 
cooling loop required.) 

2. A gas UFg tore with a liquid blanket. (External loop used to cool 
blanket . ) 

3. A gaseous UFg core with a circulating liquid blanket. 

After a brief preliminary study, the first two concepts were rejected 
arid the last one chosen for a detailed design effort. 

The first concept Was abandoned since reprocessing the blanket material 
would require a system shutdown, at least a partial reactor disassembly, 
and the installation of new blanket material. The cost of fabricating new 
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blankets would add to the total power cost. 

The second concept calls for an external coolant loop through the 
blanket which is undesirable from a neutron economy point of view. Some 
fraction of the blanket could be continuously bled off for reprocessing and 
fi$siojx product removal j a fact which makes this concept superior to concept 
1 ., 

The last concept was the one selected for the detailed design work. It 
employs a liquid blanket which will be continuously circulated for blanket 
cooling and reprocessing. No external coolant loop is required with this 
design.. The problem of finding a suitable liquid blanket material containing 
correct amounts of thorium and beryllium, while being stable and having good 
thermal-physical properties, was quickly resolved and the circulating blan- 
ket concept accepted, 

7.3 Reactors 

7.3.1 Structure and Geometry 

The initial calculations were based on a spherical core, which was 
necessary because the MACH- 1 diffusion code can only do physics calculations 
in one-dimension. A spherical reactor, however, becomes impractical due to 
fabrication and flow distribution problems. A right cylinder with ellipsoid 
heads and radius approximately equal to the blanket thickness and height 
equal to the diameter is easily fabricated and good geometry to work with 
from both a practical and a calculational point of view. 

A uniform velocity profile for the UFg flowing through the core is de- 
sirable to simplify calculations and maximvm reactor performance. Several 
methods of flattening the velocity profile were considered, but the only 
practical solution was to use many small inlets and outlets to the reactor 
core. The tubes that carry the UFg into and out of the core are offset to 
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minipize neutron losses through the tubes. 

The inlet and outlet flow passes through 100 tubes Whose total cross 
sectional area comprises approximately one tenth of the area of the end of 
the cylinder. Figure VII-2 illustrates the inlet and outlet geometry for 
the core. The tubes will constitute approximately one tenth of the blanket 
volume at the ends of the core. The introduction of UFg and the tube struc- 
tural materials into the blanket will have a degrading effect on the breed- 
ing ratio. There was no other apparent method which would distribute the 
flow aiid have less impact on the breeding ratio. 

The blanket will be pressurized to the same pressure as the core (on 
the order of 100 atm.). The core liner will thus be designed to withstand 
a pressure difference of only 10 atmospheres. The outside pressure vessel 
will need to be capable of containing the lOO atmospheres of normal operat- 
ing pressure plus a 20 % safety margin, or 120 atmospheres total. These 
pressures are not extreme and can easily be accommodated. The reactor pres- 
sure vessel is designed according to specifications from the ASME Boiler and 
Pressure Vessel Code; Section Ill-Rules for Construction of Nuclear Power 
Plant Components.^ The cylinder will have ellipsoid heads with a 2:1 mayor 
to minor axis ratio. The core liner thickness was selected as 1.27 cm. 

Since the core physics calculations indicated a height and diameter of 109.6 
cm Was necessary in order to achieve criticality, the core vessel had outside 
dimensions of 112.1 cm diameter and 130.4 cm maximum length. 

Reactor physics studies indicated that the optimum blanket thickness 
for a gas-core breeder reactor is of the order of 100 cm. Consequently, the 
inside dimensions of the reactor pressure vessel were approximately 312 cm 
diameter and 330 ciU maximum length. Because the modified Hastelloy-N pres- 
sure vessel was at the same temperature as the blanket material (1600°R), 
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the Diaxlmuin allowable stress In the pressure vessel is 10,000 p si in order 
to insure a long pressure vessel operating life (30 years). Using the tech- 
niques outlined in Reference 1, the required pressure vessel wall thickness 
is 29.3 cm. If the pressure vessel temperature is reduced by cooling, the 
wall thickness can be significantly reduced. 

The 1,27 cm wall . thickness. .for the core liner allowed 229 psi pres- 
sure difference between the UFg and blanket material to be sustained for 
long term reactor operation. In case of a rapid depressurization of the 
blanket, the core liner can withstand a pressure difference of 1,375 psi for 
a period of 6 minutes. 

All stress calculations were based upon a 30 year operating life for 
the reactor. The difference in densities between the UFg and the molten 
salt blanket produced a buoyancy force on the core liner. This force was 
calculated to be 7,660 lb g. 

The reactor and primary coolant systems will be set in a pre-strassed 
concrete reactor vessel (PGRV) . This PCRV will provide both safety and 
shielding. PGRV design was not included in this report since the details 
would depend on a precise system configuration and size. 

Both the core liner and pressure vessel should be designed to prevent, 
or at least minimize, flow stagnation areas. The core vessel uses large 
numbers of entrances and exits to minimize stagnation and the blanket ves- 
sel uses symmetry and smooth, rounded surfaces to eliminate flow stagnation. 

7.3.2 Control Systems 

In a gas-core reactor, temperature, pressure, and flow rates can be 
expected to have significant effects on core reactivity. These effects will 
be even more pronounced in this reactor than in a water-moderated, fixed- 
fuel reactor. Some type of control system must be provided to account for 


reactivity changes, as well as take care of start-ups and shut-downs. 

* 

The control system for this reactor is provided by six control drums 
evenly spaced around the core. The control drums are in the blanket * but 
as close to the core as possible to maximize their worth. This type of con- 
trol systems has been used successfully on the SNA? and NERVA reactors in 
the United States and on the TOPAZ reactor in Russia. 

Control rods of the type commonly in use ott power reactors were con- 
sidered but rejected. They would have to be inserted into the blanket instead 
of the core, and would thus have a smaller reactivity. In addition, station- 
ary guide tubes through the blanket would reduce the breeding ratio. The 
ends of the reactor design are also rather complex and many extra control 
rod penetrations would complicate the design. 

The control drums to be used on this reactor might consist of a fixed 
beryllium cylinder around which a sleeve of thermal neutron absorber (such 
as cadmium or boron carbide) rotates. By having the sleeve rotate Instead 
of the entire drum, the response time can be decreased due to the much smaller 
inertia of the sleeve. 

The control drums work by changing the size of the reflector area and 
thus the number of thermal neutrons reflected back into the core. When the 
poison sleeve is facing away from the reactor core, the beryllium cylinder 
thermalizes neutrons and reflects thermal neutrons back into the core. As 
the sleeve is moved around between the drum and the core, the neutrons strik- 
ing the poison sleeve are absorbed, effectively reducing the number of neu- 
trons reflected back into the core. 

The control drums should be designed so that any f put can shutdown the 
reactor. Under normal operating conditions, four drums will be used for re- 
actor start-ups and shut-downs and reactivity control. The other two would 


be used in aix emergency situation, such as the fai^ of a regular drum 
or during a scram., 

J.3*3 - Seals 

All of the major parts are joined together by either welding or by 
double ”0” ring seals* The "O" ring seals used here are of the same type 
most current reactors , including PWR* s , BWR* s , and the molten salt 
reactor li The *'0*- rings should be made of a metal compatible with the sub- 
stances they may come in contact with. The seal consists of two concentric 
”0’* rings, each shaped like a hollow torus j with several holes drilled around 
the inside circumference of the ring. As pressure increases, the "0" rings 
inflate, forming the seal. (See Figv VII-3) 

7*4 System layout 

Figures VII“4 and VII-5 illustrate side and top views for the proposed 
1000 Mwth power plant. The plant contains two turbines, each of which has 
a high pressure and low pressure component operating on the same shaft. 

Two regenerators are used at the exit of the low pressure turbine in order . 
to keep the size of the components down to reasonable values* The entire 
power plant will be housed inside a prestressed concrete pressure vessel 
(PCRV) which will act as the containment structure for the power plant. 

PGRV*s are presently employed as the containment for high temperature gas- 
cooled reactofs which operate at higher temperatures (1870°R) than the UFg 
gas cpfe reactot. 

7.5 Mass in Power Plant 









A - UFg Reactor p - Regenerator G - Pump 

B - HP Turbine E - Condensor H - Emergency Turbi 

C - LP Turbine F - Reheator I - Blanket Reprocessing Plant 

FIG. VII-4 LAYOUT-SIDE VIEW 
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TABLE VII -1 

MASS OF UFg IN POWER PLANT COMPONENTS 


MASS n bm) 

REACTOR 2,200 
REHEATERS (TWO) 8;900 
REGENERATORS (FOUR) 250.500 
CONDENSERS (TWO) 19,T00 
TURBINES (TWO H-P, TWO L-P)* 


PIPING* 


TOTAL — — — 

280,700 


WERE NOT 'CALCULATED 


k 


loi 


The reheaters consisted of shell-and-tube heat exchangers with counter- 
flow of the molten salt in the tubes and gaseous in the shell side. 

The tubes were of 0.5 in,. l.D. and 0^625 in. O.D. The UFg velocities were 
limited so the pressure loss across the shell side did not exceed the 0.5 
bar pressure loss assumed in the thermodynamic calculations. 

The regenerators also consisted of shell-and-tube heat exchangers with 
cpunterflpw of the high pressure liquid and boiling UFg in the tube side 
and low pressure gaseous UFg in tha shell side. Since most of the UPg mass 
was due to the liquid phase in the tubes, the tube diameters weie kapt small 
with 0.125 in. I, D, and 0.25 in* 0*D. UPg velocities were also limited so 
that pressure losses did, not exceed 0.5 bar in the shell side pass. 

The condensers Were of shell-and-tube construction with cross flow of 
the condensing UPg over the tubes. Tube diameters were 0.5 in. I.D* and 
0*625 in, O.D* 

The heat transfer calculations were performed using the standard 
techniques for analyzing heat exchangers discussed in Reference 2. In the 
regenerator, substantial temperature differences occurred in the UFg so the 
analysis was performed over these regions so that temperature variations 
in thermal -physical properties could be accounted for. 

As shown in Table VII-I, the mass of UFg in the core and heat exchangers 
Was 280,700 Ibn,, This is an unacceptably high value. The reason for the 
large mass is that the low thermal conductivity of UFg caused small heat- 
transfer coefficients Which in turn required large heat transfer surface 
areas in the heat exchangers in order to meet the required heat load. Sig- 
nificant reductions in UFg mass cOuld be achieved by putting small amounts 
of a high conductivity gas such as helium in the UFg. In addition, no 
effort was made to minimize the UFg mass in the power plant. By permitting 


higher UF^ pressure losses in the heat exchangers, the flow velocities can 
he ihcreased with subsequent increase in heat-transfer coefficients. Another 
factor is that employing heat transfer enhancement techniques, such as using 
finned tubes i used in the design of compact heat exchangers can also produce 
significant reductions in heat exchanger volumes with consequent reduced UFg 
masses. By simultaneousiy ■using all of these techniques to minimize the IlEg 
mass it is likely that an order of magnitude reduction in system mass can be 
achieved. 

Areas of Further Study • 

Many aspects of this reactor power plant require further study. Many 
problem areas became apparent only as the design progressed. Much research 
and experimental work remains before these designs can be finalized. 

The heat transfer calculations were performed using standard empirical 
correlations obtained from experiments with water and gases such as helium 
and air. Experiments need to be performed measuring heat-transfer coeffi- 
cients of UFg for flow of single-phase liquid, boiling UFg, and gaseous UFg 
in order to confirm the validity of using eonventlonal heat-transfer coeffi- 
cent correlations . 

The control drums are an area that will require a great deal of further 
study. The required sizes and materials, as Well as the corresponding re- 
activity worths, needs to be computed. The effects of the drums on the breed- 
ing ratio is also unknown. It is possible that these drums will not possess 
enough negative reactivity, considering their position out of the core, to 
control the reactor. They may also have too long a response time, moving 

I 

from full out to full in. If for either of the above two reason^ control 
drums cannot be used, then a new system will need to be devised. 
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The control drums will need to be cooled. This can be accomplished 
by drilling vertical holes In the beryllium cylinder and pumping heliuai 
through the lower control drum drive assembly into the drum. It would 
then flow upward through the holes and the gap around the edge and . out the 
control drum drive assembly at the top of the drum. The heat generation 
rates in the drums needs to be calculated, to see if helium can indeed do 
the cooling and at what rate it will need to be pumped. 

: Since the blanket is a molten salt with a high melting point, some type 

of electric heaters will need to be installed in the blanket to keep the salt 
fluid during prolonged shut-downs. The pipes in the blanket circulation 
System will also need to be electrically heated to prevent freezing. The 
feasibility needs to be analyzed , 

A thermal shield will be necessary between the outer pressure vessel and 
the PCRV. It will probably need to ba halium-cooled and heat generation rates 
in this shield need to be analyzed. 

Thermal stresses in the pressure vessel and core liner need to be analyzed 
for this system. Piping and vessel thicknesses need to be calculated. The 
core vessel and control drums may need more support than currently specified. 
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Table Al . Thermodynamic Properties of UFg as an Overheated Vapor 

specific voltime, Vi cu ft/lbjjj^ 
enthalpy, h, Btu/lbjj^ 

entropy, S, Btu/lbjjjr °S. 


T(“R) 

V 

h 

S 

V 

h 

■ s 

V 

h 

s 


p 

= 14.504 

psl 

p 

= 72.52 

psi 

p 

= 145.04 

psi 











720 


44.075 

0.3127 


44.505 

0.3042 




900 

1.828 

60.845 

0.3329 

0.386 

60.63 

0.3236 

0.172 

60.2 

0.3190 

1080 

2.235 

78.690 

0.3514 

0.441 

78.26 

0,3418 

0.215 

71.83 

0.3373 

1260 

2.634 

97.610 

0.3675 

0.522 

96.965 

0.3583 

0.257 

96.325 

0.3535 

1440 

2.882 

115.67 

0.3813 

0.600 

115.369 

0.3724 

0.298 

114.595 

0.3678 

1620 

3.390 

135.235 

0.3934 

0.675 

134.59 

0.3845 

0.337 

134.165 

0.3800 

1800 

3.760 

153.51 

0.4038 

0.752 

153.51 

0.3951 

0.375 

153.08 

0.3905 

1980 

4.151 

172.43 

0,414 

0.830 

172.43 

0.4049 

0 .414 

172.215 

0.4005 

2160 

4.528 

191.565 

0.4238 

0.905 

191.479 

0.4142 

0.453 

191.178 

0.4099 

2340 

4.902 

210.7 

0.4328 

0.982 

210.485 

0.4228 

0.491 

210.055 

0.4184 

2520 

5.285 

229.62 

0.4408 

1.058 

229.491 

0.4309 

0.530 

229,405 

0.4266 

2700 

5.660 

249.185 

0.4486 

1.134 

249.185 

0.4387 

0.535 

249.056 

0.4343 


P 

«= 290.08 

psi 

P 

=435.12 psi 

P = 

= 580*16 psi 


720 


tm M- 


s mm 

MS 

mm. 

MM 

MM 

MM 

900 

.0758 

58.695 

0.3144 

0.0418 

56.545 

0.3106 

.0232 

55.04 

.3080 

1080 

,1030 

76.97 

0.3328 

.0651 

76.11 

0.3303 

.0461 

74.82 

.3279 

1260 

.1253 

95.675 

0.3494 

.0810 

95.245 

0.3472 

.0597 

94.6 

.3453 

1440 

.1469 

113.95 

0.3638 

.0966 

113.305 

0.3610 

.0715 

112.875 

.3590 

1620 

.1673 

133.515 

0.3762 

.1108 

132.87 

0.3735 

.0827 

132.526 

.3717 

1800 

-1871 

152.65 

0.3868 

.1242 

152.306 

0.3845 

.0929 

152.005 

.3827 

1980 

.2072 

171.785 

0.3968 

.1380 

171.441 

0.3944 

.1035 

171.14 

.3927 

2160 

.2267 

190.92 

0.4061 

.1512 

190.576 

0.4037 

.1135 

190.318 

.4022 

2340 

.2459 

209.84 

0.4144 

.1644 

209.41 

0.4123 

.1236 

209.195 

.4106 

2520 

.2657 

229.19 

0.4226 

.1783 

228.975 

0.4204 

.1336 

228.803 

.4189 

2700 

.2844 

248.927 

0,4302 

.1900 

248.755 

0.4281 

.1431 

248.669 

.4264 


P 

= 725.2 psi 

P 

= 1160.32 

psl 

P = 

1740.48 

psi 

900 

.0081 

52.245 

.3045 

. 0064 



.0060 

MM 

■MM 

1080 

.0393 

73.1 

.3255 

.0177 

69.23 

.3201 

.0102 

64.93 

,3152 

1260 

.0466 

94.17 

.3439 

.0271 

92.02 

,3397 

.0169 

88.795 

.3348 

1440 

.0565 

112.445 

.3576 

.0342 

111.37 

.3543 

.0222 

109.865 

.3512 

1620 

.0658 

132.225 

.3703 

.0411 

131.15 

.3670 

.0272 

129.86 

,3641 

1800 

.0742 

151.575 

.3814 

.0447 

150.801 

.3782 

.0308 

150.07 

.3754 

1980 

.0827 

170.925 

.3914 

.0518 

170.28 

.3885 

.0346 

169.85 

.3860 

2160 

.0910 

190.146 

.4007 

.0570 

189.544 

.3977 

.0385 

188.985 

.3952 

2340 

*0990 

208.98 

.4091 

.0623 

208.55 

,4062 

.0422 

208.12 

.4038 

2520 

,1070 

228.674 

.4174 

.0674 

228.33 

.4145 

.0457 

228.115 

.4121 

2700 

.1146 

248.54 

.4350 

.0725 

248.325 

.4224 

.0490 

245.53 

.4200 
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t(°R) V- h S 

P = 2320.64 psi 


V 


V 


h 


S 


L! 



Q 

e 

D 


h 


P « 2900,8 psi P = 4351.2 psi 


900 

.0058 



.0057 

38.055 

.2858 

.0053 

38.27 

.2839 

1080 

.0102 

61.705 

.3113 

.0074 

60.845 

.3094 

.0065 

60.2 

.3070 

1260 

.0124 

86.86 

.3320 

.0103 

85.57 

.3298 

.0080 

83.85 

.3266 

1440 

.0165 

108.79 

.3571 

.0134 

107.5 

.3464 

.0098 

105.35 

.3427 

1620 

.0200 

128.785 

.3614 

V0162 

128.94 

.3598 

.0116 

126.42 

.3563 

1800 

.0233 

149.425 

.3734 

.0189 

148.651 

.3714 

.0133 

147.705 

.3686 

1980 

.0263 

169.205 

.3838 

.0214 

163.775 

,3820 

.0150 

167.7 

.3791 

2160 

.0292 

188.426 

.3932 

.0238 

187.91 

.3915 

.0167 

187.265 

.3887 

2340 

.0319 

207.905 

.4020 

.0260 

207.561 

.4005 

.0183 

206.83 

.3977 

2520 

.0348 

227.9 

.4103 

.0283 

227.556 

.4087 

.0199 

227.04 

.4061 

2700 

.0374 

247.893 

.4180 

,0805 

247.514 

.4164 

.0216 

247.25 

.4139 




Table A2. Thermodynamic Properties o£ Saturated 

D 

specific volume of saturated liquid, v^, cu ft/lbjj^ 
specific volume of saturated vapor, v * cu ft/lbj^j 
enthalpy of saturated liquid, h^, Btuflb^j^ 
enthalpy of saturated vapor, h , Btu/lbj^ 
entropy of saturated liquid, s|, Btu/lbjj^- ”R 
entropy of saturated vapor, S , Btu/lbj^-“R 


T(“R) 

P (psi) 


Sf 


h 

g 

S 

S 


^fg 

606.888 

21.756 

0 

.2389 

.004378 

35.26 

.2970 

.5595 

35.26 

626.688 

30.458 

2.494 

.2428 

.004480 

36.894 

.2978 

.4712 

34.40 

644.688 

42.352 

4.73 

. 2461 

.004566 

38.614 

.2987 

.4148 

33.88 

662.688 

55.115 

7.052 

.2499 

,004657 

40.076 

.2997 

.3420 

33.02 

680.688 

72.520 

9.374 

.2529 

.004750 

41.624 

.3003 

.2618 

32.25 

698.688 

88.474 

11.696 

.2564 

.004852 

42.871 

.3010 

.2195 

31.18 

716.688 

111.681 

14.018 

.2597 

.004959 

44.118 

.3017 

.1734 

30.10 

734.688 

137.788 

16.34 

.2625 

.005088 

45.58 

.3023 

.1394 

29.24 

752.688 

166.796 

18.92 

.2656 

.005218 

46.87 

.3028 

.1145 

27.95 

770.688 

203.056 

21.50 

.2693 

.005368 

48.16 

.3039 

.0929 

26.66 

788.688 

239.316 

24.08 

.2720 

.005531 

49.45 

.3042 

.0806 

25.37 

806.688 

282.828 

26.66 

.2755 

.005745 

50.525 

.3051 

.0683 

23.86 

824-688 

333.592 

29.412 

.2786 

.006067 

51.557 

.3055 

.0592 

22.14 

842.688 

389.432 

32.164 

.2821 

.006197 

52.374 

.3061 

.0461 

20.21 

860.688 

451.074 

35.088 

.2855 

.006669 

52.718 

.3060 

.0375 

17.63 

878.688 

523.594 

38.184 

.2889 

.007198 

52.374 

.3051 

.0284 

14.19 

896.688 

609.168 

41.624 

.2938 

.008283 

50.654 

.3038 

.0204 

9.03 

905.688 

667.184 

46.612 

.2987 

.011496 

46.612 

.2987 

.0115 

0 
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appendix b 

THERMAL-PHYSICAL PROPERTIES OEUEg 

Uranium hexafluoride is a clear, colorless, crystalline solid that 
sublimes at atmospheric pressure directly to the vapor state; Because of 
its extensive use in gaseous diffusion enrichment plants and fuel reprocess- 
ing plants, its physical and chemical properties at low and moderate temper- 
atures have been extensively studied. Some UFg properties listed in the 
( 12 ) 

literature' * ^ are given in Table Bl, 


Table Bl. Some Properties of Uranium Hexafluoride 


Property Value 


triple point, at 1134 mm Hg,“C 64.052 

sublimation point, °C 56.4 

density 

solid, g/cm^ 5.09 

liquid, g/ml 6.63 

heat of formation, solid, at 25‘’C, kcal/mole -516 
heat of vaporization, at 64.01 “C, Icoal/mole 6.907 

heat of fusion, at 64.01 “C, kcal/mole 4.588 

heat of sublimation, at 64.01 °C, kcal/mole 11.495 

critical temperature, °C 230^2 

critical pressure, atm 45.5 


Uranium hexafluoride is a highly reactive substance, acting as a moder- 
ately powerful fluOrinating agent. It reacts vigorously with water to form 
uranyl fluoride (UO 2 F 2 ) aud hydrogen fluoride. It does not normally react 
with oxygen, nitrogen, or dry air. UFg is rapidly reduced by hydrogen gas. 

UFg gas does not start disassociating until temperatures of about 1500”K 
(2700“F) are reached. Figures Bl and B2 illustrate the calculated composition 
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Eig. B“2. Composition of UFg Species at 10 Atmospheres Pressure 
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of UFg speqies at temperatures to 4000 “K for pressures of one and ten 
atmospheres, respectively. Figure B3 illustrates the ratiu of specific 
heats out to 3000 ®K, and Figure B4 shows the results of Soviet calculations 
of the composition of UFg species out to 51,000°K. Gas turbine power plants 
will probably operate only at terapetatures below 1500 °K because of corrosion 
and other problems at higher temperatures. 

Data on the ratio of specific heats is necessary for turbine and com.- 
pressor performance calculations. These data were measured to 1500“K at 
the University of Florida, and the results are presented in Figure B5. 

Viscosity of liquid UFg between 70 and 210*’C is given by the following 

( 6 ) 

aquatxon ; 


n 1 fi-7 .. (554-0. 023P)/T 

T] - 1.67 X 10 e 


(Bl) 


where P is in pounds per square inch, T is in degrees Kelvin, and T) is in 
poise. 

The viscosity of UFg vapor has been investigated by Llewellyn; and 
Fowler. Cohen summarized their data using the following empirical formulae; 


For Fowler’s data: 

Tl = 2.46 X 10“^ (85 to 165“C) 

11 1,78 X lO"® (150 to 260“G) 

For Llewellyn’s data; 

11 « 2.10 X lO"^ vp>779 200 “O 

where T is in degrees Kelvin; 1] is in poise. 


(B2) 

(B3) 


(B4) 
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RATIO OF SPECIFIC HEATS 
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Fig. B-5. Ratio of Specific Heats for Uranium Hexafluoride 






The most reliable thermal conductivity'" values of UFg vapor for 
temperatures between 0 and 100 "C are represented by the following equation ; 

= 1.46 (1 + O.OQ42t“C) X lO""^ cal/cm-sec-“C (B5) 

where is the average of the experimental and calculated conductivity 
value. 
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APPENDIX C 

THERMO -PHYSICAL PROPERTIES OP MOLTEN salt reactor PtEL 

Molten salt reactor fuel is composed of LiP (71.7 mole%), 56^2 
mole%), ThP^ (12 mole%) and UP^ (0.3 mole%). Some of its properties are^^^: 
molecular weight = 6A 
melting point = 930°F 
vapor pressure at 1150°P < 0.1 mm Hg 
Cp ** 0.324 Btu/lb. 

The fuel density and viscosity can be represented as functions of 
temperature by the following relations ; 

p(lb/£t^) = 235.0 - 0.023l7t (“F) (Cl) 

jt(lb/ft-hr) = 0.2637 exp (7362/T° (R)) (C2) 

Thermal conductivities at three different temperatures are listed below: 

at 1300 °F K = 0,69 Bt'a/lir-ft-'’F 

at 1175 °F K - 0.71 

at 1050“F K = 0.69 
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D.l Mathematical Model 


The energy equation in cylindrical coordinates is given by ; 


^ pO^T = -(7.pCpl3) - <V.g) - (T!7U) + (^) ^ 


DC 

'ot 


+ pT + q'*' 


Certain assumptions can be made to simplify this expression: 

(1) assume steady state conditions! 

(2) neglect viscous forces; 

(3) assume velocity in r and 9 directions negligible; 

(4) assxmie constant pressure in the core; 

(5) assume axial conduction negligible. 

Noting Vq s= - V*(KgVT)» Eq. (D,l) reduces to 


V.pC TCt = V.K„VT + q' 

P Li 


Expanding Bq. (D.2) and using the assumptions above 


bz p 


(pO Tu^(r,z» “ I ^ (r Kj, If) + q'" (t.z) 


and rearranging 


K, 


^ bz r br \ pC br/ ^ pC„ 




(D.l) 


(D.2) 


(D.4) 


(D.5) 




Fi 


Ihe expression for the eddy diffusivity for heat transfer is assumed 





t " (?■ r (^ - (^) ) 


Then the right side of Eq. (D.5) can he expanded 




r ^ ('<«■%> b?) “ -1J+ ^ ■ 2 (r^ )} S 


Thus Eq. (D.5) becomes 


u 


z bz 




Equation (D.8) is the mathematical model \diose solution gives the tem- 
perature profile of the design reactor. The partial derivatives may be 
estimated using finite differences 


^ ^ = ^i+i.r-^i-1.1 . ^ ^ ^i+i.i~*-"iq^^i-ia . g. 

bz S ’ br 2h ’ jjj-2 1 j2 


,-T. 


2„ T. ,-2T. ,+T. 


In finite difference form and rearranging Eq. (D.8) becomes 


T. , + ^ • 

.,2 L,1 u L 




h2 


■] 


•T. .-T. 




z o 


. S , q '» 

pC u 
p z 


(D.IO) 


where 


S = mesh size in the axial direction 
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h = mesh size in the radial ditect ion 

« = k/pCp 

A = (iJ,/p)(:72K| f)t 



Rewriting Eq* (D*10) using the heat generation term developed 


S(offGjj) 



h2 -1 

O '* 


J + rl L 2h J 


SQ 


lit 




(D.U) 


where 


u 

z 



therefore, given the core geometry, inlet temperature, power level, and 
flow rate the temperature profile of the core can be determined. 


D.2 Computer Model 

Equation (D.il) was solved using a inarching technique down the channel 
on the Georgia Tech GDC CYBER 74 computer. The program steps the tempera- 
ture radially and down the channel. Two cases of boundary conditions 
governing the solution were considered; 

Case I: (a) Uniform inlet temperature distribution (T^ = 1036. 77°R) 

(b) Insulated boundary ^ — 0 along tlie entire length of 

the core 
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W ii< Ti |l w 


bT 


(c) Sj^etrical distribution radially r— 
entire length of the core 0^ 


r«0 


= 0 along the 


Case II; (a) Uniform inlet temperature distribution - 1036, 77“R) 

(b) Insulated boundary up to the axial position 'where » 

1660 ®R; constant wall temperature thereafter 

<W = 

(c) Syiimetrical diStributioh radially. 


Equation (D. 11) poses problems when u^ = 0 (whe^ r-r^) or when r “ 0. 
Therefore, wall and centerline temperatures must be computed using another 
method. The method Used is finding an equation for an appropriate curve 
which will satisfy the boundary conditions and extrapdlate to the desired 
temperature. 

At the Wall the terapetature is assumed of the form 



i 2: 3 

*0 + V + “z" 


CD.12) 


Using the boundary condition 
points hear the wall 


bT 



0 and knowing the temperature at two 


where 


br 


= 0 - 2aj^rQ + 3a2rQ , f = r 


Ti = a^ + aj_r^^ + a2ir/ , r r^ 


^2 “ ^0 ^ 1^2 *^ 2^2 * ^ 


r^ = radius at the wall 
r^,r 2 = two points neat the wall 
Tj^,T 2 = temperatures at r. and r 


<D.13) 


Ihe unknown, coefficients axe fouud using determinants. Knowing them, T(r=XQ) 
can be readily determined using Eq. (D.I2). 

At the channel centerline the tenq»erature is assumed of the form 


T = Sq + aj_r 


(D.14) 


bT 


irtiich satisfies the boundary condition ^ 


0 . 


r=0 


Using known temperatures at two other points near the centerline 


= ao + aj^r^ , r = r^ 

*^2 “ ^0 ®1^2 ’ ^ “ ^^2 (D.15) 

The unknown coefficients are found using determinants. Then the centerline 
temperature at r = 0 follows from Eq. (D.14), i.e., T = Sq. 

D.3 Iteration, Convergence, and Stability 

An iterative approach was taken so that properties were computed near 
the (j4l) temperature as opposed to the properties being evaluated at the 
(j) temperatures. Only one iteration was needed because temperatures did 
not change significantly for any given row. 

Convergence was the most time consuming portion of the numerical analy- 
sis. Due to the complexity of Eq. (D.ll>, the effect of various parameters 
was somewhat masked, and trial -and-error was almost the rule in determining 
mesh size for convergence. Next to the stability criteria (to be discussed 
next), the size of the step in the r direction is the most significant fac- 
tor in converging to a solution. Ideally the minimum Ar should be determined 
below which the temperature profile is no longer dependent on the size of 
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Ar (withiu an acceptance criterion). The step in the r direction was 
reduced to .16 cm. It is not known if this is the minimum step size since 
smaller steps result in either exceeding the time limit on the GDC CYBER 74 
authorized for this group or, if too small, results in an instability the 
cause of which is unknown. 


Az 


Az 


The stability criterion was determined to be ^ < I ; there ^ J 


was chosen for all calculations. 
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